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i»,  AonR.cT  ’Numerical  procedures,  are  developed  for  the  digitaL solution  of  the  integral  equations  ! 
for  the  current  induced  on  a  perfectly  conducting,  two-dimensional  periodic  surface  of  arbitrary  | 
profile  when  a  plane  electromagnetic  wave  is  incident.  By  using  Floquet's  theorem  the  range  of  i 
integration  is  reduced  to  a  single  period,  and  special  summation  techniques  consisting  of  a  Poisson 
summation  and  the  subtraction  of  the  dc  term  are  used  to  improve  the  convergence  of  the  infinite 
series  representation  of  the  Green's  function.  The  integral  equations  are  then  solved  numerically  j 
using  the  moment  method  and  an  interpolation  scheme.  1 

Data  are  obtained  for  both  the  surface  ^nd  far  fields  for  a  variety  of  sinusoidal,  full-  ! 
wave  rectified,  inverted  full-wave  rectified  and  triangular  profiles  for  plane  waves  of  either  polar-! 
ization  at  oblique  as  well  as  normal  incidence,  and  the  results  are  compared  with  the  predictions  I 
of  physical  optics.  - - 

The  numerical  results  are  used  to  illustrate  some  interesting  physical  phenomena, 
notably  the  P-type  and  S-type  Wood  anomalies  associated  with  the  frequency  and  angular  responses 
of  diffraction  gratings,  and  to  develop  a  scheme  to  estimate  back  scattering  from  a  sinusoidal 
surface  at  oblique  incidence. 

The  knowledge  gained  in  the  study  of  scattering  from  periodic  surfaces  is  taen  applied 
to  the  study  of  rough  surfaces  by  treating  the  surface  as  a  small  scale  roughness  superimposed 
upon  a  periodic  base  (representing  the  large  scale  roughness).  The  small  scale  roughness  is 
approximated  by  a  random  function  with  a  Gaussian  distribution. 
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Chapter  I 
INTRODUCTION 


Periodic  structures  find  a  variety  of  applications  in  fields  rangjing  from 
acoustics  and  electromagnetics  to  optics.  Typical  examples  are  the  diffraction 
gratings  used  in  optical  spectrometers  and  as  resonator  components  in  sub- 
optical  lasers,  the  slow-wave  structures  in  microwave  tubes,  and  for  micro- 
wave  lenses  and  broadband  absorbers.  They  are  also  important  to  oceanograph¬ 
ers  in  the  study  of  acoustic  or  low-frequency  radio  wave  scattering  from  an 
idealized  sea. 

Mathematically,  the  problem  of  the  scattering  of  waves  by  a  periodic 
surface  is  one  which  is  easy  to  formulate  but  difficult  to  solve.  It  consists 
of  solving  the  wave  equation  subject  to  boundary  conditions  imposed  at  a  sur¬ 
face  whose  shape  ca.n  be  very  complicated.  For  this  reason,  an  exact  analytical 
treatment  of  the  problem  is  seldom  possible,  and  it  is  only  recently,  with  the 
advent  of  the  high  speed  computer,  that  reliable  results  have  been  obtained. 

Prior  to  this  time,  a  variety  of  approximate  analj'tlcal  treatments  were 
developed,  most  of  them  based  on  an  approach  originated  by  Lord  Rayleigh 
(1878).  Assuming  the  surface  to  be  Infinite  in  extent  and  making  use  of  the 
periodicity  in  (say)  the  coordinate  x,  Rayleigh  expanded  the  scattered  field  in 
a  discrete  spectrum  of  outgoing  plane  waves,  which  representation  was  assumed 
to  hold  right  down  to  the  surface.  Application  of  the  boundary  condition  leads 
to  an  infinite  set  of  linear  algebraic  equations,  valid  for  all  x,  from  which  to 
determine  the  complex  amplitudes  of  the  scattered  waves,  infinite  in  number. 

The  crucial  and  unjustified  step  (see  Lippmann,  1953)  in  this  procedure 
is  the  assumption  that  the  representation  of  the  scattered  field  holds  every¬ 
where  above  and  on  the  boundary.  This  has  ever  since  been  a  subject  of  con¬ 
troversy  (see,  for  example,  Lysanov,  1958).  Thus,  for  a  perfectly  conducting 
sinusoidal  surface,  Rayleigh  (1878)  obtained  a  solution  by  successive  approxi¬ 
mation  based  on  the  initial  neglect  of  all  attenuated  waves  and  assuming  that 
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the  amplitudes  of  the  corrugations  are  small  irregularities  compared  to  the 
wavelength.  Later,  for  the  small  surface,  Tai  (1948)  proceeded  via  an 
orthogonal  mode  expansion,  followed  by  matrix  truncation,  and  others  have 
pursued  essentially  the  same  path.  An  analogous  treatment  for  a  corrugated 
interface  between  two  homogeneous  media  was  developed  by  Rayleigh  (1907),  and 
applied  to  a  sinusoidal  profile  by  Pavageau  (1963)  and  to  a  triangular  profile  by 
Bousquet  (1963). 

A  refinement  of  Rayleigh's  method  was  developed  by  Uretsky  (1965), 

His  method  starts  in  the  samvs  way  as  Rayleigh’s,  but  instead  of  merely  assuming 
that  the  scattered  field  can  be  expanded  as  an  infinite  set  of  plane  waves,  Uretsky 
proved  that  this  is  possible  for  observation  points  not  too  close  to  the  boundary. 
To  avoid  Rayleigh’s  assumption  that  the  representation  holds  right  down  to  the 
surface,  Uretsky  used  Green’s  theorem  to  express  the  scattered  field  as  an 
integral  over  the  elementary  sources  induced  on  the  boundary  by  the  incident 
wave.  Numerical  results  were  obtained  for  a  sinusoidal  surface  with  a  Dirichlet 
boundary  condition. 

Methods  of  small  perturbations  similar  to  Rayleigh's  were  used  by 
Miles  (1954)  and  Katsenelenbaum  (1955).  The  typical  procedure  is  as  follows: 
the  boundary  conditions  specified  on  the  uneven  surface  y  *  f{x)  are  trans¬ 
formed  to  the  plane  y  =  0  by  expansion  in  a  power  series  with  respect  to  f(x). 
Thus  the  problem  of  scattering  from  an  uneven  periodic  surface  is  I’educed  to 
solving  the  wave  equation  subject  to  specified  boundary  conditions  on  a  plane. 

At  best,  all  such  solutions  are  valid  only  for  corrugations  whose  height 
is  much  smaller  than  the  free  space  wavelength,  and  in  an  attempt  to  overcome 
this  restriction,Meecham  (1956a)  used  a  variational  method  to  find  the  angular 
distribution  of  scattered  energy  for  a  perfectly  conducting  grating.  The  scat¬ 
tered  field  was  represented  as  a  linear  combination  of  known  solutions  of  the 
wave  equation  whose  coefficients  were  obtained  by  a  least  square  fit  to  the 
boundary  condition  at  the  surface,  and  the  procedure  was  then  applied  to  a  tri¬ 
angular  (or  sawtooth)  profile. 
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A  somewhat  different  approach  was  tftteen  by  Eckart  (1853),  He  expressed 
the  scattered  field  aa  a  Helmholts  integral  involving  the  scattered  field  itself  and 
it«  normal  derivative  on  the  boiodary.  The  integral  was  then  evalasted  after 
assuming  that  the  reflection  took  place  at  an  equivalent  fiat  surface^  and  approxi*' 
mating  the  scattered  field  and  its  normal  derivative  on  the  boundary  by  the 
incident  field  al«me.  Another  me^od  also  based  on  variatitxui  on  the  Kircfahoff 
approximatieQ  was  used  by  Brekhovskildi  (19S2).  The  method  assumed  that  the 
field  at  the  uneven  surface  was  local  in  nature  and  could  be  specified  in  terms 
of  the  laws  of  geometrical  optics^  The  field  at  each  point  on  the  surface  was 
assumed  to  be  the  same  as  if  the  reflection  from  that  point  were  to  occur  from 
an  Infinite  plane  tangent  to  the  surface  at  the  npecified  point. 

A  method  based  on  irisyslcal  optics  was  used  by  Senior  (1959).  Once 
again  the  determination  of  the  field  is  reduced  to  quadratures,  and  Senior  Aowed 
that  for  a  plane  wave  at  normal  incidence  on  a  sinusoidal  grating,  the  physical 
optics  integral  can  be  evaluated  exactly  to  give  the  complex  amplitudes  of  the 
scattered  waves.  As  indicated  in  Appendix  A,  the  same  is  true  (if  shadowing  Is 
■  Ignored)  for  both  polarizations,  and  for  oblique  incidence  as  well  as  normal.  It 
should  be  emjdiasized,  however,  that  the  solution  is  still  approximate  by  virtue 
of  the  poetulated  surface  field  distribution,  and  the  failures  of  the  physical  optics 
estimate  of  the  surface  fiolds  are  examined  in  Chapters  m  and  V  of  this  work. 

A  method  which  is  quite  distinct  from  all  of  the  above  was  developed 
by  Slvbv  (1964)  who  used  conformal  transformation  and  a  consideration  of  the 
static  limit  to  analyze  reflections  from  periodic  surfaces  with  shallow  and  deep 
corrugations.  The  procedure  is  similar  to  that  recently  employed  by  Millar 
(1969,  1971)  to  investigate  the  inherent  limitation  of  Rayleigh's  method  (see 
Chapter  vn).  As  first  noted  by  Llppmann  (1953),  it  is  not  in  general  valid  to 
assume  that  the  expansion  of  the  scattered  field  as  a  discrete  spectrum  of  out¬ 
going  waves  alone  holds  over  the  entire  scattering  surface,  and  this  fact  was 
later  verified  by  Petit  and  Cadilhac  (1964)  in  the  case  of  a  sinusoidal  grating. 
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In  any  general  treatment  of  the  grating  problem  it  is  therefore  necessary  to 
allow  "ii^oing"  waves  in  the  immediate  vicinity  of  the  surface. 

All  of  the  methods  considered  so  far  are  very  limited  (satisfactory  only 
for  certain  classes  of  structures)  and  have  little  success  in  predicting  certain 
characteristics  of  periodic  structui’es  (for  example  the  P-type  Wood  anomalies). 
In  the  desire  for  a  better  understanding  of  Wood  anomalies  (see  Chapters  IV  and 
VI),  a  different  point  of  view  was  proposed  by  Hessel  and  Oliver  (1965).  They 
used  a  guided  wave  approach  and  replaced  the  uneven  surface  by  a  pMne  which 
had  a  modulated  surface  reactance.  Their  theoiry  is  more  general  than  the  pre¬ 
vious  theories  and  predicts  the  general  location  and  shapes  of  the  anomalies. 
Since  the  surface  geometry  and  the  modulated  plane  reactance  cannot  be  related 
quantitatively,  information  about  the  performance  of  a  specific  structure  cannot 
be  obtained. 

Recently,  another  approach  based  on  function  theoretic  techniques  was 
used  by  Tseng,  Hessel  and  Oliner  (1968),  and  Itoh  and  Mittra  (1969).  Unfor¬ 
tunately  structures  for  which  such  methods  are  possible  are  very  few. 

Without  exception,  all  of  these  analytical  attempts  to  determine  the  fields 
scattered  by  periodic  surfaces  are  subject  to  approximation,  either  Implicit  or 
explicit,  and  it  is  only  with  the  use  of  high  speed  computers  permitting  the 
direct  digital  solution  of  the  integral  equations  that  reliable  results  have  been 
obtained. 

Most  of  the  initial  work  in  this  area  was  carried  out  by  the  French 
investigators,  such  as  Petit,  Cadilhac  and  Wirgln,  and  was  motivated  by  the 
desire  for  more  efficient  optical  diffraction  gratings.  In  his  early  papers  Petit 
(e.g.  1963)  followed  Rayleigh's  approach  In  expressing  the  scattered  field  as  a 
discrete  spectrum  of  outgoing  waves  alone,  leading  to  a  matrix  equation  for  the 
determination  of  the  spectral  amplitudes.  Since  the  matrix  was  then  truncated 
and  inverted  numerically,  it  will  be  appreciated  that  the  method  is  no  more  than 
a  digitization  of  that  originated  by  Rayleigh.  However,  in  later  papers  (Petit, 


1967),  the  Rayleigh  assumption  was  circumvented  by  using  integral  equation 
formulation.  Series  expansions  were  adopted  for  the  tacidant  and  scattered 
fields  and  the  Integral  equation  converted  to  a  matrix  equation  which  was  solved 
numerically.  Specific  results  were  obtained  for  plane  wave  incidence  on  gratings 
with  triangular  profiles,  and  the  efficiencies  computed.  A  rather  dlffereiit 
approach  was  taken  by  Pavageau  (1967)  who  dferived  the  integral  equation  direct¬ 
ly  in  terms  of  the  unknown  surface  current.  The  equation  waa  cast  in  the  form 
of  a  hCHihomogeneous  Fredholm  equation  of  the  second  kind  and  solved  by  iteration. 

Methods  which  are  very  similar  to  that  which  we  shall  use  have  recently 
been  employed  by  Neureuther  and  Zaki  (1968,  197ia,b)  andliy  Oreen  (1970).  The 
former  considered  scattering  by  periodic  structures,  either  dielectric  or  per¬ 
fectly  conducting.  The  integral  equation  was  obtained  from  Green’s  theorem 
and  the  modified  Green's  function  expressed  in  either  of  two  ways  depending  on 
the  parameters  of  the  surface.  The  first  (space  harmonic)  representation  is 
analogous  to  that  used  by  Petit;  the  second  consists  of  an  infinite  series  of  Haakel 
functions,  and  was  computed  using  a  Mellin  transform  and  an  asymptotic  com¬ 
parison  scheme.  Numerical  data  were  obtained  for  sinusoidal  profiles,  and  very 
recently  data  were  obtained  also  for  triangular  proffjvt'A  by  Kalhor  and  Neureuther 
(1971).  Green  (1970)  also  used  the  space  harmonic  r'opresentation,  but  improved 
its  convergence  by  summation  tedmlques.  Data  war*'  presented  for  the  surface 
field  and  diffraction  efficiencies  of  perfectly  coiKluctkig  gratings  with  triangular 
profiles. 

The  present  work  also  employs  the  numerical  approach,  and  Is  concerned 
with  the  scattering  of  electromagnetic  waves  by  Infinite;  perfectly  ooiiducting, 
two-dimensional  periodic  surfaces  of  arbltraiy  but  coi-'toious  profile.  Plane 
wave  incidence  is  assumed,  with  either  E  or  H  polarlvnua-.i,  and  both  noxmal 
and  oblique  incidence  are  considered.  A  representation  of  the  .Jreen's  function 
is  employed  which  is  sLnilar  to  that  used  by  Green  (1970),  iiiid  the  convergence 
is  improved  still  further  by  subtracting  the  dc  terms.  Tbia  has  the  added  ad¬ 
vantage  of  making  explicit  the  behavior  of  the  Green's  function  In  a  neighborhood 


of  its  singularity.  The  integral  equations  are  formulated  on  the  basis  of  the 
Green’s  theorem,  and  the  range  of  integration  is  reduced  to  a  single  period 
by  application  of  Floquet’s  theorem.  Using  the  moment  method  and  an  inter¬ 
polation  scheme,  the  integral  equations  are  converted  into  a  matrix  equation 
which  is  solved  to  yield  the  unknown  irarface  current. 

Chapter  n  is  concerned  with  the  formulation  of  the  integral  equations 
and  the  de-  ''’•'pment  of  numerical  procedures  for  solving  them.  Numerical  data 
for  E  polarization  are  presented  in  Chapters  HI  and  IV.  Chapter  m  is  concerned 
with  the  surface  field  while  Chapter  IV  is  devoted  to  the  scattered  field.  Numer¬ 
ical  results  are  obtained  for  a  variety  of  surface  profiles  such  as  sinusoidal, 
full-wave  rectified,  inverted  full-wave  rectified  and  triangular,  having  various 
combinations  of  amplitudes,  periods  and  for  various  Incidence  angles.  In  pre¬ 
senting  the  numerical  data  in  Chapter  IV,  particular  attention  is  given  to  the 
Wood  anomalies  (P-type),  and  in  both  chapters,  the  data  are  compared  with  the 
physical  optics  predictions. 

The  numerical  results  for  H  polarization  are  presented  in  Chapters  V  and 
VI  and  the  S-type  Wood  anomalies  discussed.  The  results  support  Palmer’s 
experimental  observation  that  the  S-type  anomalies  are  generally  stronger 
than  the  P-type  and  can  occur  on  structures  with  relatively  shaLow  grooves. 

Chapter  vn  is  devoted  to  Rayleigh’s  method.  After  a  brief  historical 
survey  of  the  assumption,  numerical  data  are  obtained  for  the  scattered  field 
using  the  exact  representation  as  well  as  Rayleigh's. 

Lastly,  In  Chapter  VIH,  the  knowledge  gained  from  the  study  of  periodic 
surfaces  is  applied  to  an  investigation  of  scattering  from  rough  surfaces.  The 
rough  surface  is  treated  as  a  small  scale  roughness  superimposed  on  a  periodic 
base  (large  scale  roughness)  with  the  small  scale  roughness  represented  by  a 
random  variable  having  a  Gaussian  distribution. 


Chapter  Q 

FORMULATION  OF  THE  INTEGRAL  EQUATION 
FOR  NUMERICAL  SOLUTION 


ft** 


We  oonalder  here  an  infinite,  perfectly  conducting  periodic  surface 
illuminated  by  a  plane  electromagnetic  wave.  Since  the  surface  is  assumed 
two  dimensional  in  the  sense  of  being  independent  of  a  Cartesian  coordinate  z, 
the  entire  problem  is  two  dimmaional,  and  the  most  general  solution  can  be 
deduced  from  the  particular  solutions  appropriate  to  incident  pbuie  waves  having 
either  ^  or  H  in  the  z  direction,  i.e.  parallel  to  the  corrugations.  In  either 
case,  the  problem  is  essentially  scalar. 

2.1  Formulation 

It  is  convenient  to  develop  first  the  integral  equation  in  the  somewhat 

simpler  case  of  a  scattering  surface  of  finite  extei^.  Let  S  (see  Fig.  2>1)  be 

this  suzlaoe,  and  surround  it  by  another  closed  surface  Let  0(r)  ami  G(rjr*) 

be  two  scalar  functions  which  are  continuouf;,  together  with  their  first  and  second 

derivatives,  on  S  and  S  and  throughout  the  volume  V  enclosed  by  them.  Assume, 

K 

moreover,  that  0(r)  satisfies  the  homog^yieous  wave  equation 


(V^  +  k^)0(r)  •*  0 


(2.1) 


inside  and  on  the  boundaries  of  V,  wbereas 

(V^  +  k^)G(r|i*')  ^  -  At- t')  . 

Applying  Green's  theorem  to  the  volume  V,  we  obtain 

S|^G(r(r')-~  <p(T')-i}Kr')-^  G(r|r')|  dS' 


(2.2) 


(2.3) 
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9 


by  virtue  of  the  properties  of  the  delta  function.  We  now  identify  0(r)  with  a 


total  field 


where  i//\r)  is  an  incident  plane  wave  originated  by  a  source  at  infinity  (and 
therefore  outside  V),  and  4/  (r)  is  the  field  scattered  by  the  surface  S.  Since 

d 

ip  (r)  must  satisfy  a  radiation  condition  at  infinity,  its  contribution  to  the  inte> 

gral  over  S-  decreases  to  zero  as  the  surihice  S  recedes  to  infinity,  whereas 
i  R  R 

p  (r)  contributes  itself.  Equation  (2.3)  then  becomes 

=0^(r)+  \  I'GfrjrO  g^0(r')-<^r')g^Q(rjr')  •  dS'  (2.4) 

Js 

and  in  spite  of  the  assumptions  of  an  incident  plane  wave  and  a  surface  of  finite 
extent,  Eq.  (2.4)  is  also  valid  for  an  arbitrary  incident  field  and  for  a  surface 
S  extending  to  infinity.  T.  latter  case,  howsver,  the  proof  is  by  no  means 
trivial  (Jones,  1952). 

The  particular  situation  of  concern  to  us  is  that  in  which  S  is  doubly 
infinite  and  divides  space  into  two  regions.  It  is  then  sufficient  to  integrate 
over  only  the  upper  ("illaminated")  side  of  the  surface.  We  also  assume  that 
S  is  independent  of  the  coordinate  z,  and  this  allows  uc  to  distinguish  two  par¬ 
ticular  cases  ao«>rding  to  the  polarization  of  the  Incident  field.  Mks  units  are 
emplo3red  and  a  time  factor  e"^  suppressed  throughout. 

2.1.1  E  (or  horizontal)  Polarization 

If  E^  *  the  electric  vector  in  the  scattered  field  will  also  be 
—  z 

confined  to  the  z  direction,  and  we  can  make  the  identification 


4At)  =  E*%) 
% 


(2.5) 


Since  the  normal  derivative  of  E^  is  related  to  the  surface  current  density 
by  the  equation 
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^  e‘V)  =-)>«' • 

where  n  is  the  permeability  of  free  space,  Eq.  (2.4)  can  be  written  as 


(2.6) 


\  ]  G(£li:')V£'>+ j^E*”'(r')^G(rJr')^  dS'  . 

)s  "  (2.7) 


At  the  perfectly  conducting  surface  S  the  boundary  condition  is 


E“(rt  =  0 
2  — 


(2.8) 


and  hence,  on  allowing  r  to  approach  the  surfece,  Eq.  (2.7)  gives  rise  to  the 
integral  equation 


G<r|r')K^(r')dS’  =  j^E‘(r) 


(2.9) 


For  a  two  dimensional  problem  the  free  space  Green's  function  is 

G(rlr’)  =  Hjj^\klr-r'|)  (2.10) 

(2) 

where  is  the  zero  order  Hankel  function  of  the  second  kind.  The  final 
integral  equation  is  therefore 


Vr')H®(k|r-r'|)d<' 


(2.11) 


from  which  K  has  to  be  determined. 

2 

2.1.2  H  (or  vertical)  Polarizatio.i 

If  2,  the  total  magnetic  field  is  likewise  in  the  z  direction,  and 

Xt 

on  making  the  identification 
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^r)  »  H‘'nr) . 

Eq.  (2.4)  becomes 


(2.12) 


H^(r)  =  H^(r)+' 
z  —  z  — 


G(r|r')|dS'.  (2.13) 


sJS 

The  boundary  condition  at  the  perfectly  conducting  surface  S  is 

^Hf(r)=0  (2.U) 

and  on  allowing  r  to  approach  the  surface,  Maue  (1949)  has  shown  that  (2. 13) 
reduces  to 


1 

2 


Hf(r) 


(2.15) 


The  quantity  H  is  the  induced  current  density,  aud  since  this  current  flows 
z 

tangential  to  the  surface  in  a  plane  perpendicular  to  the  z  direction,  we  can 
write 

IC(r)  *  H^°^r)t  (2. 16) 

— t-  z  - 

where  t  is  a  unit  vector  tangential  to  S.  Inserting  finally  the  expression  (2. 10) 
for  the  free  space  Green’s  function,  we  have 


i.e. 


HQ^^(k(r~r’|)df'  »  4j|^H^(r)- 1  K^(r)  j  , 


(2.17) 


This  is  an  integral  equation  from  which  to  determine  K^, 


2.2  Reduction  of  the  Integral  Equations  for  the  Periodic  Surfaces 

We  now  make  use  of  the  fact  that  the  incident  field  is  a  plane  wave  and 
that  the  surface  y  =  f(x)  is  periodic  with  period  d,  i.e. 

f(x  +  md)  =  f(x)  ,  m  =  1 1,  i’2,  1 3, . . .  .  (2.19) 

As  before,  it  is  convenient  to  consider  separately  the  cases  of  E  and  H  polari¬ 
zation. 

2.2.1  E  Polarization 
Let  us  assume  that 

(2.20) 


where  0  is  the  angle  of  incidence  with  respect  to  the  normal  to  the  mean  sur¬ 
face.  The  integral  equation  is  as  shown  in  Eq.  (2.11),  but  since  the  integration 
extends  from  -oo  to  oo,  this  is  not  appropriate  for  "  numerical  solution.  How¬ 
ever,  by  invoking  the  periodic  property  of  the  surface,  the  integral  can  be  reduced 
to  one  over  a  single  period  alone  at  the  expense  of  a  more  complicated  form  for 
the  kernel. 

From  Floquet's  theorem,  we  have  that 


K  (r  +  mdx)  =  K  (r)e 

7."  7,— 


-jkmdsin© 


m 


=  ^1,  t2.  +3, 


(2.21) 


Moreover, 

d?'  =  +  dx' 

where  f'(x)  is  the  derivative  of  f  with  respect  to  x.  This  allows  us  to  express 
the  integral  as  one  along  the  x'  axis,  and  since 
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ll-r'l  *  \/^x-x*-md)^+ (y-y')^  , 


we  obtain 


f  J]]  (k{(x  -  X*  -  md)^  +  (y-  y')" 

Jo 


2  ^  -jkmdsinO 


/l+  {f(x')} 


2 


4  -Jk(x8ln©-ycos0) 
"  0 
uyi 


K  (x')dx' 
z 


(2.22) 


valid  for  0  <  X  <  d  with  y  “  f(x),  y’  »  f(x'). 

In  arriving  at  <2.22)  we  have  implicitly  assumed  that  interchanging  the 
order  of  integration  and  summation  is  valid  (a  fact  which  is  no  means  obvious), 
but  even  so  Eq.  (2.21)  is  not  a  very  promising  equation  for  numerical  purposes 
because  of  the  extremely  ppor  convergence  of  the  Infinite  series  constituting  the 
kernel.  To  rectify  this  situation,  consider 


x‘-md)^+(y-y’)^ 


m=*oo 


-jkmdsine 


(2.23) 


Using  the  Fourier  integral  representation  of  the  Hankel  function,  the  Poisson 
summation  formula  (see,  for  example,  Morse  and  Feshbaoh,  1953)  applied  to 
(2.23)  gives 


2  1  -J( -~+k8in0)(x-x')  -j|y-y'|x^ 

Pi-3  2_.  3r“ 

m=-(X)  m 


(2.24) 


where  _ _ _ 


(2.25) 


and  the  chosen  branch  of  the  square  root  is  that  for  which 

ImX  <0 
m  " 


In  order  to  satisfy  the  radiation  condition. 


The  expression  (2.23)  for  has  been  used  by  Yen  (1962)  and  Green 
(1970),  and  is  equivalent  to  that  obtained  by  Neureuther  and  Zaki  (1968)  using 
a  space  harmonic  expansion.  The  features  that  should  be  noted  are  that  if 
y  ^  y'  the  successive  terms  decrease  exponentially  with  increasing  (m,  for  all 


'•  a.  I  o 

m+  —  sin0 
2ir 


kd 

27r 


-Iv 


and  that  in  addition  there  is  an  algebraic  decrease  (proportional  to  m  )  pro¬ 
vided  by  the  factor  X  in  the  denominator. 

m 

We  can  produce  still  a  further  improvement  in  the  convergence  properties 
of  the  series  by  separating  out  the  zero-frequency  (k  =  0)  terms,  and  this  has 
the  added  advantage  of  making  explicit  the  behavior  of  in  the  vicinity  of  the 
singular  point  x'  =  x,  y'  =  y.  Putting  k  =  0  in  (2.24),  we  have 


00 


-j-^(x-x')  -•^|m(y-y')| 


2it 


P  ' 

‘lk=o 


2 


d  1 1  I  27r 

m=-(X)  j|mj  -r 


(2.26) 


and  since 


00 


m=l 


t- 


m=-oo 


-  -i— (x-x')-~iv-v' 


,  27r ,  ,v  2tt  .  I 


(see,  for  example,  Collin,  1960),  provided  y'  f  y,  it  follows  that 


15 


P  =  (m*0  term)- ^  log  <1-6 
^  k*0  ^  U- 


-f  Oy-y’l-J<*-*']i 


(2.27) 


Hence,  aubtnustlng  the  zexo  frequency  value  of  each  term  of  (2.24)  and  then 
adding  (2.27),  the  integral  equation  (2.21)  takea  the  form 

[  Gjte, y;  X', y)K^(x')el'“'*‘“®  {f(x-)}='  ax’  - ^  , 

Jo 


0<x<d , 


(2.28) 


where 


G  (,  .  _i_  ,-lk|y-y'looBe 

'»j>x,y,x,y)  a 


Id,  Tr  r  -^Qy-yl-Kx-x'^TI 


-J^(x-x-)  -^ly-yr 

d  ^  e  .  e 


+  7  .  e 


^  f -)|y-y1x„  -^ly-y'l'] 


m“l 


(2.29) 


»  yk^(‘^~  +  ksind)^  ,  X~  “  \fk^-  ("~^-  keln0)^  .  (2.30) 
m  u  m  '  d 


I 


If) 


We  observe  that  the  Green's  function  has  now  been  expressed  as  the  sum  of 
a  term  involving  explicitly  the  amplitude  and  the  angle  of  incidence,  a  logar¬ 
ithmic  function  representing  the  true  behavior  of  G  in  a  neighborhood  of  its 
singularity,  and  two  series  which  are  themselves  differences  of  two  convergent 
series. 

2.2.2  H  Polarization 

The  incident  field  in  this  case  is  taken  to  be 


„i-  ,  -jk(x  sine -y  cos  e) 

H  (r)  =  e 
z  — 


(2.31) 


anH  the  reduction  of  the  integral  equation  (2. 17)  to  a  form  suitable  for  numerical 
solution  proceeds  in  much  the  same  way  as  for  E  polarization.  Since 

a  a  I 


a 

an' 


- 1 


f  ax'  "  ay'  /  ’ 


yi  +  -(_f'{x'); 

Eq.  (2.17)  can  be  written  as 

.d 

P2K^(x')dx'  =  4j  -^H^r)  -  ^  K^(r)  j 

Jo 

where 

00 

P  «  X!  H®  (k-,(x-x'-md)^+(y-y’)^ 

m=-oo 


(2.32) 


'  ^-jkmdsine 


00 


-  ^  y  -L 

d  4-J  X 


,  f'(x')(^+ksine)-X^sgn(y~y')  x 


m=-oo  m 


A  e 


.j(^+ksine)(x-x')  -j|y-y"|X^^ 

e 


(2.33) 


in  which  sgn(T)  is  the  signum  function  and  is  given  in  (2.25).  The  zero 
frequency  Ik  =  0)  limit  of  is 


:,3^isc2ss@rg^'^-; 
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00 


=  1  IZ|-^J'(x')+i6gn(y-y')[ 

ik  -0  m=-oo  ‘  j 


.  2mT[  ,  ..  Ztt  ,  . 

-3  — {x-x')  -“|m(y-y'), 

e  e 


00  ^  2x0.11  ,  ,  I 

2,'  ^  9  -  "3 — -x')-  y-vM 

=  -  ^f'(x‘)+jsgn(y-y')^  {f'(x')+3  sgn(y-y’)|  ^ 

m=l 


00  .  2x011  r 

m=l 


and  since 


00 


r. 

m=l 


2  '  2  ^^^**2  * 


we  have 


\  -  |f'<x')  +  ^-[f'(x')+j8gn(y-y')|  cotj^-  J  (x-x ')+j  ~  |y-y' ]  j 

+  ^  j-f'(x')+jsgn(y-y')|cot|^  (x-x')  +  j|iy-y'l^  .  (2.' 


34) 


Hence,  by  subtracting  the  zero-frequency  value  of  each  term  of  (2.33)  and  then 
adding  (2.34),  the  integral  equation  (2.32)  takes  the  form 

G2(x,y;x',y')K^(x')dx' =Jd|2H^(x)-K^(x)-  ,  0<x<d,  (2.35) 


where 
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G_(x,y;x',y’) 


=  -j  sgn(y-y')- j -'|^f'(x')tan0-sgn(y-y')|  e  jocose  |y  y 


00  .  2miT . 

_  .  _,-jksme(x  -  X')  V  ■ 


-jly-y'(x^ 


X 


m 


2m7r 


|y-y'l 


^2ni7r  +  H 

'  f'(x')  — j-  +  kam0)-x^sgii(y-y')> -e  {jf'(x')- sgn(y-y')j- 


^  .2m£  r  -j|y-y'|X 

.^-jksine(x-x')^g^  d  ^  ■- 

m=l  X 


m 


i_  m 


•  f'(x')  +ksin0)-X  sgn(y-y’)^  + 

V  \  <1  m  ' 


2rmr  . 

ly-y'i 

+  e  i  if'(x' 


jf'(x')  +  sgn(y-y')C 

J 


+  '|-|  f'(x')+jsgn(y-y')j-  cot  ^  (x- x')+j  [y- y' 

+  2  ■^3 ■^ly-y'i  > 


(2.36) 


In  spite  of  the  obvious  complexity  of  the  Green's  function  G^,  its  form 
is  directly  analogous  to  G^^  in  consisting  of  certain  explicit  terms,  a  pair  of 
cotangent  functions  representing  the  true  behavior  of  Gg  in  a  neighborhood  of 
its  singularity,  and  two  series  which  are  themselves  differences  of  two  con¬ 
vergent  series.  Once  again,  therefore,  the  singularity  of  the  Green's  function 
has  been  separated  out. 

We  further  note  that  the  integral  equations  (2.28)  and  (2.35)  have  been 
derived  without  approximation.  In  consequence,  the  formulation  so  far  is 
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exact  for  any  two  dimensional  perfectly  conducting  periodic  surface  which  is 
smooth  in  the  sense  of  having  a  continuous  first  derivative. 

2 .  ■;  Comments  on  the  Green*8  Functions  of  the  Integra!  Kquations 

Before  we  discuss  the  numerical  procedures  of  solving  the  integral 
equations  shown  in  Eqs.  (2,28)  and  (2.35),  it  is  appropriate  to  examine  their 
kernels  (modified  Green's  functions)  with  some  care.  Equation  (2.28)  for  E 
polarization  is  a  Fredholm  Integral  equation  of  the  first  kind  with  a  weakly 
singular  kernel.  Very  little  is  known  about  equations  of  this  type,  and  even  the 
uniqueness  of  the  solution  is  not  assured.  On  the  other  hand,  Eq.  (2.35)  for  H 
polarization  is  a  Fredholm  integral  equation  of  the  second  kind  with  a  weakly 
singular  kernel.  The  advantage  of  this  equation  over  that  for  E  polarization  is 
that  we  can  use  an  iterative  scheme  (Neumann  series)  to  obtain  the  solution. 

This  has  been  done  by  Payageau  and  Bousquet  (1969)  using  the  geom.6tric  optics 
expression  for  the  surface  field  as  a  first  approximation.  Ad  shown  in  Appendix 
B,  however,  the  Integral  equation  for  E  polarization  can  be  converted  from  one 
of  the  first  kind  to  one  of  the  second  kind,  so  the  successive  approximation 
method  can  be  applied  here  as  well. 

If  we  examine  the  Green's  function  for  E  polarization,  Eq.  (2.29),  it  is 
found  that  the  surface  sloptJ  dependence  of  the  kernel  occurs  explicitly  only  in 
the  terms  \'l+f’(x)  .  For  surfaces  of  gentle  slope  such  that  |f'(x)|  «  I,  it  is 
now  obvious  that  the  surface  current  will  be  relatively  insensitive  to  slope .  The 
second  term  of  the  Green's  function  has  a  logarithmic  singularity  which  repre¬ 
sents  the  correct  singular  behavior  for  E  polarization  in  two  dimensions.  The 

last  two  terms  consist  of  infinite  series.  Whenever  or  X  Is  purely 

m  m 

Imaginary,  the  corresponding  terms  are  exponentially  damped.  Hence,  in  a 
numerical  solution  of  the  integral  equation,  these  series  can  be  terminated  at 
a  term  somewhat  beyond  that  where  the  exponential  decay  starts.  On  physical 
grounds,  it  would  appear  that  as  functions  of  K,  the  values  of  these  series  could 
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change  abruptly  when  an  individual  mode  becomes  cut-off.  This  could  then 
lead  to  an  abrupt  change  in  surface  field  distribution.  However,  for  surfaces 
ol  small  periods  (d/X  <  0.2),  the  contribution  of  the  infinite  series  is  insig¬ 
nificant,  and  the  first  two  terms  dominate. 

For  H  polarization,  the  Green's  function  is  much  more  complicated  and 
more  dependent  on  angle  of  incidence  and  surface  slopes,  but  its  qualitative  be¬ 
havior  is  not  too  different  from  that  for  E  polarization.  If  f'(x)  =  0,  the  surface 
is  a  flat  sheet  and  it  can  be  shown  that  vanishes .  The  current  distribution 

is  now  given  by  the  explicit  term  alone,  and  this  is  just  the  geometrical  optics 
value.  Again  the  correct  singular  behavior  for  H  polarization  in  two  dimensions 
is  represented  by  the  last  two  cotangent  functions. 

2.4  Numerical  Procedures 

Methods  for  the  numerical  solution  of  integral  equations  have  been 
extensively  discussed  in  the  literature  (see,  for  example,  Harrington,  1968). 

The  general  procedure  consists  of  reducing  the  equation  to  a  finite  set  of  alge¬ 
braic  equations,  i.e.  to  a  matrix  equation,  and  can  be  illustrated  by  considering 

C' 

\  G(x,x')K(x')dx'=  F(x)  ,  0<x<d.  (2.37) 

Jo 


We  assume  that  the  unknown  function  K(x')  can  be  expanded  in  terms  of  linearly 
independent  base  functions  0^(x')  such  that 


K(x') 


a  0  (x') 


(2.38) 


where  the  a 
gives 


n 


are  the  associated  constants. 


Substitution  of  (2.38)  into  (2.37) 


d 

G(x,  x')0^(x')dx'  =  F(x) 

0 


(2.39) 
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and  the  solution  of  the  integral  equation  has  now  been  reduced  to  the  determination 

of  the  constants,  a  ,  n  =  1, 2, . . N. 
n 

There  are  several  possible  ways  of  finding  the  a  ,  e.g.  least  squares  fit, 

n 

Galerkin's  method,  and  the  collocation  method,  and  it  is  the  last  of  these  that 
we  shall  employ. 

The  collocation  method  converts  (2.39)  into  a  system  of  N  linear 
equations  by  forcing  the  two  sides  of  (2.39)  to  be  equal  at  N  sampling  points  in 
the  interval  (0,  d).  This  is  simply  a  point- matching  procedure  and  results  in  the 
matrix  equation 

N  pd 

Ea  \  G(x  ,x')0  (x')dx'  =  F(x  ),  0<x  <  d,  m  =  1,  2, . .,  N  ,  (2.40) 

n\  m  n  m  —  m-  >  »  > 

Jo 


There  now  remains  the  problem  of  choosing  the  base  functions  i5^(x'),  and  here 
again  there  are  several  possible  choices,  e.g.  ractangular,  quadratic  and 
sinusoidal.  By  appropriate  choice,  we  can  economize  in  the  number  of  sampling 
points  requred  for  an  "accurate"  approximation  to  the  solution  K(x),  and  exper¬ 
ience  has  shown  that  a  rectangular  function  is  not  in  general  a  good  choice, 
whereas  sinusoidal  interpolation  often  works  rather  well. 

The  particular  form  of  sinusoidal  Interpolation  that  we  have  adopted  is 
predicated  on  the  use  of  sampling  points  which  are  uniformly  distributed  in 
0  <  x'  <  d.  The  range  of  integration  is  therefore  broken  up  into  N  increments 

of  length  A  'N  (see  Fig.  2-2).  Furthermore,  let  x'  be  the  midpoint  of  the 

1 

n’th  cell,  i.e,  x'  =  (n-  -  )A,  and  let  Ax'  denote  the  interval 
~  i  n 


x'  ■ 
n 


<x'  <  x'  ■ 
2  -  -  n 


n 

_A 

2 


We  assume  that 


A  +B  sink(x'-x  )  +  C  cosk{x'-x  ),  ifx'eAx' 
n  n  n  n  n  n 


K(x')  = 


\ 


(2.41) 


0 


otherwise 
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and  specify  the  constants  A^,  and  by  continuing  the  form  appropriate  to 
x'  e  Ax'  out  to  the  centers  of  the  adjacent  cells  and  imposing  continuify.  Defining 


K  =  K(x  )  and  K  .  =  K(x  ,  , we  have 

n  n  nil  nil 


K  =A  +C  , 
n  n  n 


K  .  =A  +B  sinJtA+C  coskA  , 
n+1  n  n  n 


K  =A  -B  sin,cA-hC  coskA, 
n- 1  n  n  n 


from  which  we  obtain 


A  = 
n 


-K  +2K  coskA-K  , 
n+1  n _ n^l 

2(coskA-  1) 


B  = 


K  -  K  , 
n+1  n-1 


n  2sinkA 

K  -2K  +K  , 

P  _  n+1  n  n-1 

n  2(coskA-l) 

Substitution  of  (2.41)  and  (2.42)  into  (2.40)  now  gives 

r  rcos  kA  -  cos  k(x' -  x_) 

\  G(x  ,x')  K  - r-r — r- 

I  m  n  coskA- 1 

“=‘Jax  -  - 


n 


(2.42) 


+  K 


n+l 


■-sinkA+(coskA- l)sink(x' -X  )  +  sin kA cos k(x' - x  ) 

_ _ S - S_ 

2  BinkA(coBkA- 1) 


-sinkA-(co8kA- l)sink(x' -X  )  +  sinkAcosk(x* -x  ) 


+  K 


n-1 


2 sin kA( coskA  -1) 


dx' 


=  F(x  )  , 
m 


(2.43) 
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where  x  =  (m  -  ~  )A,  m  =  1, 2, . . N.  We  note  in  passing  that  when  n  =  1 
m  £■ 

or  N  the  periodicity  of  the  problem  must  be  used  to  determine  the  constants 

K  ,  or  K  required  in  (2.43). 

-1  irt-1 

The  above  procedure  is  immediately  applicable  to  the  integral  equations 

for  E  and  H  polarizations  on  inserting  the  appropriate  values  for  the  Green's 

function  and  the  forcing  function^  but  a  brief  comment  is  desirable  concerning 

the  treatment  of  the  singular  cell.  The  Green's  functions  of  concern  to  us  are 

singular  when  x  =  x',  the  singularity  being  logarithmic  for  E  polarization,  and 

a  first  order  pole  for  H  polarization,  and  it  is  therefore  necessary  to  modify 

the  numerical  scheme  when  x  =  x  .  In  line  with  the  usual  practice,  we  divide 

m  n 

the  singular  cell  into  three  portions: 


,  A  e  . 
(x  X  , 

n  z  n  2 


(x  -  T  ,  X  +1 )  . 
n  4  n  2 


Vf' vl> 


with  0<  €  <  A.  The  first  and  third  segments  are  handled  by  the  standard 
numerical  technique,  whereas  the  central  portion  is  treated  analytically  by 
means  of  a  limiting  process  (Andreasen,  1964). 

It  must  be  mentioned  here  that  the  sinusoidal  interpolation  by  no  means 
issumes  that  the  surface  current  distribution  is  sinusoidal.  In  feet  other  methods 
based  on  quadratic  or  flat-top  functions  can  be  used  just  as  well.  However, 
numerical  results  indicate  that  good  convergence  is  obtained  by  using  sinusoidal 
interpolation  in  this  case. 

In  order  to  make  the  computer  program  applicable  to  periodic  surfaces 
of  arbitrary  shape,  it  is  convenient  to  represent  the  surface  by  its  Fourier  series. 

If  y  f(x)  is  a  periodic  function  of  period  d,  then 


“o 


00 


i: 

n=i 


^2mi  .  .  f  2m  \  ' 

— X,  xj^ 


a  cos 
n 


(2.44) 
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where 


It  is  known  from  the  theorems  of  Fourier  series  that  if  f(x)  is  a  periodic, 
continuous  and  piecewise  continuously  differentiable  function,  then  the  Fourier  series 
(2.44)  will  converge  ’.miforinly  to  fix).  Subject  to  these  restrictions  on  f(x), 
we  can  therefore  synthesize  an  arbitrary  periodic  surface  to  any  desired  degree 
of  accuracy.  In  general,  more  Fourier  terms  are  needed  for  H  polarization 
than  for  E  polarization  to  obtain  the  same  degree  of  accuracy  because  the  Green's 
function  in  the  former  case  involves  more  terms  with  the  first  derivative. 

2.5  Convergence  and  Accuracy  of  the  Numerical  Solutions 

The  numerical  solution  of  the  integral  equations  for  both  polarizations 
has  been  programmed  for  an  IBM/67  computer  in  Fortran  IV  language.  In  solving 
an  integral  equation  numerically  it  is  always  desirable  to  carry  out  a  convergence 
test  for  the  following  reason?* 

a)  to  determine  the  number  of  sampling  points  necessary  to  achieve 
the  required  accuracy,  and 

b)  to  test  whether  or  not  the  numerical  solution  approaches  a  stable 
value  as  the  number  of  sampling  points  is  increased. 

Although  it  has  not  been  shown  mathematically  that  the  accuracy  can  be 
improved  by  increasing  the  number  of  sampling  points,  it  would  appear  reason¬ 
able  to  believe  so  on  physical  grounds.  In  general  it  is  considered  that  a 
numerical  solution  is  "satisfactory"  if  it  remains  essentially  unchanged  when 
additional  terms  or  modes  are  included  in  the  computation,  but  there  are  ex¬ 
ceptions  to  this  (see,  for  example,  Lee,  Jones  and  Campbell,  1971). 
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The  accuracy  of  the  numerical  solution  in  this  work  has  been  checked 
using  the  following  criteria: 

a)  Energy  conservation.  This  furnishes  a  common  check  on  the 
accuracy  of  the  solution  in  electromagnetic  scattering  from  perfectly 
conducting  surfaces.  This  is  essentially  the  check  employed  by  Petit 
and  Cadilhac  (1964),  Neureuther  and  Zaki  (1968)  and  Green  (1970)  as 
well  as  many  others.  However,  as  pointed  out  by  Amitay  and  Galindo 
(1968)>  energy  conservation  does  not  provide  a  measure  of  accuracy  of 
a  solution  found  by  the  Ritz  or  other  related  methods,  and  hence,  in 
order  to  use  energy  conservation  to  check  other  than  computational 
round-off  errors,  care  must  be  exercised  in  choosing  the  method  of 
solving  the  integral  equation  numerically. 

The  method  that  has  been  adopted  in  this  work  is  a  subsection  and 
sinusoidal  interpolation  scheme,  and  since  the  solution  does  not  auto¬ 
matically  satisfy  energy  conseivation,  we  are  able  to  use  this  as  a 
check.  The  relation  that  must  be  satisfied  is 


^Re(X  )=kcos0 
m 


(2.45) 


where  N  is  the  number  of  propagating  modes  and  the  are  the 
scattered  field  amplitudes.  These  are  derived  in  Chapters  IV  and  VI 
for  E  and  H  polarizations  respectively. 


b)  Convergence  test.  The  numerical  solution  is  found  for  the  same 
surfece  using  different  numbers  of  sampling  points.  As  shown  later, 
the  numerical  solution  approaches  a  stable  value  as  the  number  of 
sampling  points  increases,  which  value  can  be  considered  as  the  solution. 

It  Is  found  that  the  number  of  sampling  points  necessai-y  to  yield  accurate 
results  depend  on  many  factors,  such  as  the  angle  of  Incidence,  the  ampli¬ 
tude  and  period  of  the  profile,  the  smoothness  of  the  profile  as  well  as 
polarization  of  the  incident  fields.  In  general,  for  both  polarizations 
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better  convergence  is  obtained  for  normal  incidence  or  when  the 
incidence  angle  is  small.  The  computer  program  adopted  in  this  work 
is  based  on  uniform  sampling  cells,  but  it  may  be  true  that  a  nonuniform 
sampling  scheme  would  ^jive  faster  convergence  for  the  same  number  of 
cells.  Because  the  Green's  functions  are  sensitive  to  a  change  in  the 
surface  slope,  particularly  for  H  polarization,  it  is  necessary  for  us  to 
insert  more  sampling  points  into  those  regions  where  the  surface  slope 
changes  most. 

c)  Flat  surface  check.  A  very  simple  and  useful  check  on  the  numerical 
solution  is  to  let  the  surface  become  flat.  We  should  then  obtain  the  geo¬ 
metrical  optics  value  for  a  flat  sheet,  and  this  is  indeed  the  case. 

d)  Check  with  tho  result  obtained  by  other  authors.  It  is  always 
desirable  to  be  able  to  check  the  results  obtained  with  those  obtained 
by  other  workers  (and  other  methods).  The  numerical  results  obtained 
here  will  be  compared  with  those  obtained  by  Zaki  and  Neureuther,  and 
by  Green.  Good  agreement  is  indicated. 

2.6  Computational  Time 

Perhaps  one  of  the  most  serious  restrictions  imposed  on  numerical  work 
is  the  computational  time  involved.  Although  the  computer  programs  have  been  cast 
in  a  general  form  applicable  to  surfabes  of  any  size,  economies  dictate  that  only 
reasonable  sizes  (say  from  0.2\  to  2\  in  period)  are  explored.  It  is  understood 
that  most  of  the  computational  time  consumed  is  used  in  filling  up  the  matrix 
elements.  The  rest  is  used  in  the  matrix  Inversion  process.  In  general,  the 
program  for  H  polarization  takes  more  computational  time  than  for  E  polarization, 
because  of  its  more  complicated  Green's  function.  To  save  computational  time, 
separate  programs  have  been  written  for  the  special  case  of  sinusoidal  surfaces. 
However,  in  spite  of  the  fact  that  for  tne  case  of  normal  incidence  symmetry  can 
be  used  to  reduce  the  matrix  size  by  a  factor  of  two,  and  hence  to  reduce  the 
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computational  time  by  a  factor  of  four,  separate  programs  were  not  written 
to  take  advantage  of  this,  since  most  emphasis  will  be  given  to  the  more 
interesting  case  of  oblique  incidence.  Examples  of  the  computational  time 
involved  are  given  later. 


Chapter  HI 

SURFACE  FIELD  FOR  E  POLARIZATION 

In  this  chapter  surface  field  distributions  obtained  by  numerical 
solution  of  the  integral  equation  for  E  polarization  are  presented.  The  par¬ 
ticular  profiles  considered  are  sinusoidal,  full-wave  rectified,  inverted  full- 
wave  rectified,  and  triangular  (see  Fig.  3-1)  with  various  periods  and  angles  of 
incidence.  In  each  case  the  physical  optics  approximation  is  also  presented  as 
a  basis  for  comparison. 

The  physical  optics  method  is  a  widely-used  technique  for  estimating 
the  scattered  field.  Basically,  it  is  an  extension  of  geometrical  optics  in 
which  the  surface  field  is  approximated  by  its  geometrical  optics  value,  thereby 
reducing  the  determination  of  the  scattered  field  to  quadratures.  Consequently, 
the  physical  optics  method  is  a  high  frequency  approximation  technique,  and  it 
is  suitable  only  when  the  "effective"  dimensions  of  the  scatterer  are  large 
compared  with  the  wavelength,  and  there  is  no  shadowing  involved.  How  much 
shadowing  affects  the  solution  is  still  unknown,  and  a  matter  that  is  then  of  some 
debate  is  whether  shadowing  should  be  taken  into  account,  and  if  so,  how. 

According  to  the  physical  optics  method  the  current  induced  on  the  surface 
is 

K  =  2n  ah'  in  the  illuminated  region,  and  (3.1) 

=  0  in  the  geometrical  shadow  (3.2) 

where  n  is  the  outward  normal  to  the  surface,  and  for  a  two  dimensional  sur¬ 
face  this  is  given  by 
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30 


A 

n  = 


1 


{l  +  f'{x7 


[-f'(x)x+y] 


(3.3) 


Thus,  using  Eqs.  (2.20),  (3.1)  and  (3.3)  we  obtain  for  E  polarization 


(3.4) 


in  the  illuminated  region. 

By  comparison  with  the  exact  solution,  it  is  found  that  in  general  the 
physical  optics  method  fails  to  give  accurate  results.  However,  if  we  consider 
a  "smooth"  surface  with  sufficiently  large  ratio  of  period  to  surface  height, 
we  can  then  rely  on  the  physical  optics  method. 

3.1  Surfaces  of  Small  Period  (d/X  ^0.2) 

Figures  3-2  through  3-13  represent  the  current  distributions  on 
periodic  surfaces  of  four  different  profiles  (see  Fig.  3-1).  Each  surface  has 
a  period  d  =  0.2A.,  and  a  maximum  amplitude  a  =  0.1\.  The  angles  of  inci¬ 
dence  considered  are  0=0  (normal  incidence),  0  =  30°  and  0  =  60°.  In  each 
figure,  the  exact  computed  values  are  shown  as  circled  points,  and  arc  joined 
by  a  broken  line  only  to  guide  the  eye;  the  physical  optics  approximation  is 
shown  as  a  solid  line. 

Let  us  examine  the  sinusoidal  surface.  We  observe  that  for  this 
relatively  small  period  most  of  the  current  is  concentrated  in  the  vicinity  of 
the  surface  peaks,  with  the  current  almost  zero  in  the  troughs .  There  is, 
indeed,  almost  an  exponential  decrease  in  the  current  modulus  away  from  the 
peaks,  and  the  main  effect  of  increasing  the  incidence  angle  is  to  scale  the 
curves,  leaving  the  general  shape  unchanged.  In  spite  of  shadowing  due  to 
large  angles  of  incidence,  the  curve  for  the  modulus  still  retains  its  symmetry. 


Normalized  Surface  Field  for  a  Sinusoidal  Surface 
(d  =  0. 2X,  a  =  0,  IX  and  d  =  30°)  for  E  polarization 
—0 —  exact,  —  physical  optics. 


Normalized  Surface  Field  for  an  Inverted  Full-Wave 
Rectified  Surface  (d  =  0. 2X,  a  =  0,  IX  and  6  =  0)  for  E 
polarization.  — o —  exact,  —  physical  optics. 
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Phase  (degfrees)  Normalized  Current  Modulus 


Phase  (degrees)  Normalized  Current  Modulus 


-40  I 

Fig.  3-12: 


Noriiuili/.fcl  Surface  Field  for  a  friangular  Surface  (d  (),2A.. 

(  -  0.  IX,  a  0.  IX  and  0  .'io'*)  lor  E  polarization.  — o —  exact 

—  physical  optics. 
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The  phase  is  somewhat  more  sensitive  to  0,  and  whereas  the  curve  is  almost 
flat  for  0  =  0,  the  shape  changes  noticeably  with  increase  in  0. 

Since  kp  ,  =  0. 064  (p  ,  is  the  minimum  radius  of  curvature)  for 
min  min 

the  surface,  it  is  not  surprising  to  find  that  the  i^ysical  optics  approximation 
bears  no  resemblance  to  the  exact  data.  According  to  physical  optics  the 
current  witnin  the  shadow  r^ion  should  be  zero,  but  exact  results  prove  this 
is  not  so,  and  this  is  natural  because  the  current  must  be  a  continuous  function 
along  the  surface.  The  physical  optics  phase  is  also  poor  for  0  =  0,  but  agrees 
better  for  large  0,  at  least  in  an  average  sense.  In  general  it  is  true  that  the 
physical  optics  approximation  shows  best  agreement  with  the  exact  data  in 
regions  where  the  radius  of  curvature  is  largest.  That  is,  in  all  cases  the 
maximum  discrepancy  occurs  at  the  surface  peaks  and  troughs  where  the  radius 
of  curvature  is  a  minimum. 

Examination  of  the  remaining  profiles  show  that  the  shape  of  the  modulus 
curve  is  very  similar,  except  in  small  details,  to  the  corresponding  one  for  the 
sinusoidal  surface,  and  the  phase  still  remains  relatively  constant  for  0  =  0. 

In  general,  the  above  remarks  for  the  sinusoidal  surface  apply  also  to  these 
profiles,  but  each  profile  has,  because  of  its  own.  geometry,  some  particular 
characteristics.  For  example,  in  the  case  of  the  inverted  full-wave  rectified 
profile  with  0  =  60°  (Fig.  3-7),  we  can  see  a  small  standing  wave  pattern  in  the 
phase  Just  outside  the  shadow.  This  is  probably  due  to  multiple  scattering  of  the 
incident  field  by  the  surface.  Another  example  is  the  "edge  effect"  which  is 
clearly  evident  in  the  results  for  the  triangular  profile. 

In  order  to  ^btain  a  better  understanding  of  the  surface  field,  let  us  go 
back  to  the  Integral  equation.  If  we  look  at  Eq.  (2. 29)  carefully,  we  observe 
that  for  a  sufficiently  small  period  we  can  approximate  the  kernel  by  the  first 
two  terms,  and  of  these  two  terms  only  the  first  term,  which  has  no  singularity, 
contains  the  incidence  angle.  Thus  the  approximate  integral  equation  becomes 
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(3.5) 


where 


G*(x  rx'  y')  =  — i—  - 

r  ’y*  'y’  kcoss 


-f  log|{l-e 


It- 


(3.6) 


Since  the  integral  equation  is  singular,  a  rough  estimate  of  the  solution 
to  the  integral  equation  may  be  obtained  by  examining  the  equation  in  the  neigh¬ 
borhood  of  the  singularity.  Thus 


|y-y«j— >0  iy_y'|__>0 


(3.7) 


where 


|{l-e  X 


lim  log(f  j  =  lim  log  |l-e  ^  ^ 

|x-x'|— >0  |x-x'i— >0  ^  ow  r  -1 

|y-y')— >0  |y-y'i— >0  r  - [|y-y’l  - j(x-x')J  )“| 


A  -SI  -e 


r  --j|y-y’l  --jly-y'l 

lim  log-|l+e  -2e  cos 

|x-x’j — >0  ^ 

iy-y’l— >0 

r  -y|y-y’!  --^ly-y’i  , 

lim  log -<1+6  -2e  v 

|y-y’j->o  -  J 
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Usiiig  the  following  expression. 


X 

e 


=  1 


2 

+X+I-  + 


3 

X 

31 


+. 


f 


we  obtain 


lim  log[^  3=  lim  loge  , 
Ix-x'l— >0  e— >0 

jy-y'l— >0 

where 


Thus  in  the  neighborhood  of  singularity,  Eq.  (3.5)  becomes 


(3.8) 


\^k  cos  B 


(3.9) 


where  |X-x'j<6.  6>0i8an  arbitrarily  small  number.  If  we  assume  that  the 
“ikx^sin0  /  2 

factor  K  (x’)e  *  yl+f(x')  remains  constant  in  this  small  region,  we 

z 

can  solve  Eq.  (3.9)  to  get 

+^1  6K  (x)e-i'“=‘“®||l'f  PO?  =  .  (3.1( 

Lkco8  0  27rJ  z  ’  un 

It  now  follows  immediately  that  if  the  period  d  is  sufficiently  small,  the 
first  term  in  the  square  bracket  dominates,  and  so  the  modulus  of  the  surface 
field  is  seen  to  be  proportional  to  co8  0,  a  scaling  factor.  This  is  consistent 
with  the  numerical  results. 
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The  effect  on  the  surface  field  for  a  sinusoidal  profile  due  to  change  in 
the  amplitude  as  the  period  is  kept  constant,  and  0  =  0  is  shown  in  Figs.  3-2, 

3-14  and  3-15.  It  is  observed  that  as  the  height  decreases,  so  does  the  current 
concentration  near  the  peaks;  the  phase,  on  the  other  hand,  is  much  more  nearly 
constant,  and  is  more  akin  to  that  for  a  flat  surface  than  a  sinusoidal  one, 

3.2  Surface  of  Slightly  Larger  Period  (0.2Sd/X^  0. 6) 

Figures  3-16  and  3-17  show  the  surface  field  for  a  sinusoidal  profile 
with  d  =  0.4X,  a  =  0.2\  and  0=0,  60®  respectively.  Compar'iSon  with  Figs. 

3-2  and  3-4  indicates  that  doubling  both  d  and  a  has  little  effect  on  the  modulus 
of  the  current,  but  has  a  marked  effect  on  the  phase.  For  normal  incidence  the 
phase  curve  still  has  a  small  sag  at  the  surface  trough  region.  Results  for 
d  =  0.6X,  a  =  0.3X  and  0  =  60°  are  shown  in  Figs.  3-18  and  3-19  for  the  inverted 
full-wave  rectified  and  full-wave  rectified  surfaces  respectively.  Although  the 
overall  shape  of  the  modulus  curves  bears  much  resemblance  to  the  previous 
cases  of  small  period  (Figs.  3-7  and  3-10),  there  is  evidence  (small  undulations 
along  the  curves)  of  stronger  multiple  scattering  here.  Again,  the  physical 
optics  method  fails  to  give  satisfactory  results  in  all  cases. 

3.3  Surfaces  of  Larger  Period  (d/X  ^0.6) 

Let  us  now  consider  surfaces  of  larger  period  for  which  the  requirements 

of  physical  optics  are  more  nearly  satisfied.  Referring  to  Eq.  (2.30)  it  is  seen 

that  as  d  becomes  larger,  the  number  m  which  makes  and  X  purely 

mm 

imaginary  increases  (implying  more  propagating  modes),  and  the  infinite  series 
of  the  kernel  therefore  becomes  more  important.  One  consequence  of  this  is 
that  the  surface  current  distributions  for  surfaces  of  large  period  are  much  more 
complicated  and  quite  different  from  those  for  surfaces  of  small  period . 

Figures  3-20  and  3-21  represent  two  cases  of  a  sinusoidal  surface  at 
normal  incidence  with  d  =  0.95X  and  a  =  0.25X;  d  =  1.90X  and  a  =  0.25X 
respectively.  In  the  first,  the  modulus  displays  a  strong  oscillation  having  a 
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Normalized  Surface  Field  for  a  Full-Wave  Surface  (d  =  0.6X, 
a  =  0. 3X  and  0  =  60  )  for  E  polarization,  — o —  exact, 

—  physical  optics. 
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standing  wave  ratio  of  about  6  and  a  period  of  0. 6X.  Maxima  occur  at  a 
distance  of  0. 17X  from  the  peaks  while  a  minimum  exists  at  the  trough  of  the 
surface.  The  exact  phase  is  still  characterized  by  the  familiar  almost-constant 
value.  The  physical  optics  approximation  gives  no  agreement  with  the  exact 
solution,  but  in  the  second  example  for  which  the  period  has  been  doubled  while 
keeping  the  height  constant,  the  physical  optics  approximation  is  better,  though 
there  are  some  noticeable  discrepancies.  The  maxima  of  the  modulus  still 
occur  at  the  surlEace  peaks,  but  minima  appear  where  the  surface  slope  changes 
sign,  that  is,  where  the  surface  crosses  the  x  axis.  The  minimum  radii  of  cur¬ 
vature  for  these  two  cases  are  =  0.572  and  2.30,  respectively. 

A  case  showing  better  agreement  between  the  exact  solution  and  the 

physical  optics  approximation  is  given  in  Fig.  3-22,  where  the  surface  is  a 

sinusoidal  one  with  d  *  1.6X,  a  »  O.IX  (kp  ,  =  4.07),  and  6  =  45*^.  There  is 

min 

no  obvious  explanation  for  the  apfMrent  shift  between  the  exact  and  physical 
optics  curves. 

To  examine  the  effect  on  the  surface  field  of  a  graduaJ  increase  in  the 
surface  height  when  the  period  and  angle  of  incidence  are  kept  fixed,  we  select 
the  sinusoidal  surface  with  d  =  1. 155X,  9  =  60°,  and  a  varying  from  0.05X  to 
0.7X,  The  minimum  radii  of  curvature  for  these  cases  are  as  follows: 


tlltude  a  (X) 

^min 

0.05 

4.260 

0.1 

2.130 

0.2 

1.065 

0.3 

0.710 

0.4 

0.533 

0.5 

0.426 

0.6 

0.355 
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The  numerical  results  are  illustrated  in  Figs.  3r23  through  3-30. 

It  is  seen  that  the  agreement  between  the  physical  optics  prediction  and  the 
exact  solution  deteriorates  gradually  as  the  amplitude  is  increased.  If  we 
examine  the  phase  curves  with  a  little  care  we  find  that  the  physical  optics  phase 
agrees  better  with  the  exact  phase  in  the  illuminated  region,  and  both  of  them 
are  relatively  constant  there.  In  contrast  to  the  situation  for  surfaces  of  small 
period,  the  field  is  very  sensitive  to  a  change  in  surface  height,  and  it  is  there¬ 
fore  very  difficult  to  accurately  predict.  The  undulations  in  the  modulus  curves 
are  evidence  of  multiple  scattering  which  is  particularly  noticeable  when  a  =  0. 3X. 
If  the  surface  height  is  sufficiently  large  (a  >  0. 5X),  the  surface  field  near  the 
region  of  a  surface  trough  is  relatively  small. 

Since  the  field  is  so  small  over  most  of  the  concave  portions  of  the 
surface,  it  is  of  interest  to  see  how  much  the  surface  field  as  a  whole  is 
affected  by  the  actiml  geometry  of  the  concavities.  We  examine'  this  by  replacing 
the  concave  part  (d/4  <  x  <  3d/4)  of  a  sinusoidal  surface  by  a  flat  part,  and  the 
results  are  presented  in  Fig.  3-31  for  a  =  0.5X  and  0  *  60®.  Comparison  with 
Fig.  3-28  shows  that  the  replacement  has  little  effect  on  the  surface  field  in  the 
illuminated  region,  but  does  have  a  marked  effect  tn  fite  shatfow. 

Additiaoal  data  are  presented  in  Figs.  3-32  and  3-''^  for  a  full-wave  recti¬ 
fied  profile  where  d  =  0.75X,  a  =  0.6\  and  9  -  41.8°,  and  a  triangular  profile 
where  d  =  1.75\,  0^^  =  20°,  0^  =  66°  and  0  =  12.2°.  For  the  rectified  profile 
(Fig.  3-32),  the  current  modulus  no  longer  has  a  maximum  at  the  surface  peaks 
where  kp  is  a  minimum.  In  Fig.  3-33,  the  current  modulus  shows  clearly  the 
"edge  effect"  in  the  neighborhood  of  the  edge,  and  the  physical  optics  phase 
agrees  relatively  well  with  the  exact  phase,  except  near  the  trough  region.  In 
both  cases  the  modulus  is  zero  at  the  center  of  the  surface  trough.  Fig.  3-34 
represents  the  case  of  an  inverted  full-wave  rectified  profile  with  d  =  0.95X, 
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Normalized  Surface  Field  for  an  Inverted  Full-Wave 
Rectified  Surface  (d  =  0.  95X,  a  =  0. 25X  and  0  =  0)  for  E 
polarization,  — o—  exact,  —  physical  optics. 


Fig,  3-34; 


70 


3.4  Coupling  Among  Cells 

One  aspect  which  has  not  beeii  looked  into  so  far  is  the  effect  on  the 

surface  field  of  a  gradual  increase  in  the  separation  between  cells  while  the 

other  parameters  are  kept  constant.  To  examine  this,  we  selected  a  full-wave 

rectified  surface  with  a  =  0..3X,  d  =  0.3X,  d  =  60°  (kp  =  0.0304)  and  allowed 

min 

the  distance  t  between  two  neighboring  periods  to  increase  from  0.3X  to  0.7X 
in  O.IX  steps.  The  results  are  shown  in  Figs.  3-35  through  3-39.  In  the  first 
two  figures  the  modulus  curve  has  the  shape  characteristic  of  small  periods, 
while  the  physical  optics  phase  oscillates  about  the  exact  phase.  As  I  gets  larger, 
noticeable  changes  appear  in  the  modulus  curve.  One  reason  for  this  can  be  seen 
from  Eqs.  (2,20)  and  (2.30);  as  (  increases,  the  terms  in  the  infinite  series 
corresponding  to  higher  in's  (or  higher  orders  of  propagating  modes)  become 
more  significant,  and  have  a  greater  effect  on  the  surface  cuiTent  distribution. 

As  I  passes  through  the  value  0.535X,  X  ^  changes  from  purely  imaginary  to 
real,  implying  tliat  the  m  =  -1  mode  becomes  propagating.  We  therefore  expect 
a  rather  rapid  change  in  the  surface  field  for  f  near  0.535X,  and  this  indeed 
occurs  as  can  be  seen  by  comparing  Figs,  3-37  and  i-38.  Further  increases  in 
i.  can  again  produce  large  changes  in  the  surface  field,  especially  in  the  shadow 
region  (see  Fig.  3-39),  but  it  is  worth  noting  that  thraighout  the  illuminated 
region  ‘he  modulus  is  not  very  sensitive  to  £. 

In  all  cases  considered,  the  physical  optics  predictions  bear  no 
resemblance  to  the  exact  results,  although  the  estimate  for  the  phase  is  some 
what  better  than  that  for  the  modulus.  Examination  of  the  plots  also  reveals  the 
important  fact  that  the  coupling  among  cells  is  not  very  strong  even  for  small 
periods,  and  this  leads  us  to  conjecture  that  as  regards  the  surface  field  the 
infinite  surface  can  be  approximated  a  finite  surface  consisting  of  only  a  few 
cells. 

3.5  Accuracy  Checks  and  Computational  Time 

As  mentioned  in  section  2.5,  the  numerical  results  were  checaed  by 
considering  \)  ‘uo  convergence  of  the  numerical  solution,  b)  the  limiting  approx- 
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Fig.  3-35:  Normalized  Surface  Field  for  a  Modified  Full-Wave  Rectified 
Surface  (d  =  0.  3X,  a  =  0.  3X,  d  =  60°  and  i  =  0.  3X)  for  E 
polarization,  — o —  exact,  —  physical  optics. 
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Fig.  3-;i6:  Normalized  Surface  Field  fo’  a  Modified  Full-Wave 

Rectified  Surface  (d  =  0. 3\,  a  =  0.  3X,  6  =  60°  and 
H  =  0.  4X)  for  E  polarization,  — o —  exact,  —  physical 
optics. 
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polarization,  — o —  exact,  —  physical  optics. 
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Fig.  3-38: 


Normalized  Surface  Field  for  a  Modified  Full-Wave  Rectified 
Surface  (d  =  0.  3X,  a  =  0.  3X,  0  =  60°  and  i  =  0.6X)  for  E 
polarization  — o—  exact,  —  physical  optics. 
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imation  of  a  flat  surface,  and  c)  a  comparison  with  results  obtained  by  previous 
authors.  All  these  will  be  discussed  in  this  section. 

Figure  3-40  shows  the  result  of  a  convergence  test  applied  to  the  case 
of  a  sinusoidal  profile  with  d  =  0.2X,  a  =  0.  IX  and  0=0.  It  can  be  seen  that 
the  solution  does  appear  stable,  and  that  surprisingly  good  results  are  obtained 
even  with  so  few  as  four  sampling  points.  The  limiting  approximation  of  a  flat 
surface  is  demonstrated  with  a  sinusoidal  profile  where  d  =  0,2X,  a  =  0.  05X 
and  0=0,  and  the  results  are  shown  in  Fig.  3-41. 

Note  that  in  spite  of  this  relatively  small  amplitude  the  modulus  still 
shows  significant  deviations  from  the  value  2.0  for  a  flat  surface.  This  is 
probably  due  to  the  rather  slow  convergence  of  the  infinite  series  of  Eq.  (2,29) 
when  the  surface  height  is  very  small.  However,  the  phase  agrees  exceptionally 
well  with  the  exact  phase,  0^,  for  a  flat  surface. 

Comparison  of  the  numerical  results  with  those  obtained  by  previous 
authors  is  given  in  Figs,  3-42  and  3-43.  In  Fig.  3-42,  comparison  is  made  with 
those  by  Zaki  (196J)  for  a  sinusoidal  surface  where  d  =  1.9X,  a  =  0.25X  and 
0  =  0,  and  in  Fig.  3-42  with  Green  (1970)  for  a  triangular  profile  where  d  =  1 .75X, 
0^  *  20^,  0^  *  66°  and  0=  12.  2°.  In  both  cab  s  good  agreement  is  obtained. 

The  computational  (CPU)  time  used  in  each  computation  depends  on  the 
following  factors;  a)  the  number  of  sampling  points,  b)  the  number  of  Fourier 
terms  used  to  represent  the  periodic  surface,  c)  the  number  of  terms  used  in 
U  p  infinite  series  summation  and  d)  the  accuracy  of  the  numerical  integration 
required.  In  each  case  the  accuracy  of  the  numerical  integration  using  the 
Newton-Cotes  method  was  restricted  to  1  per  cent.  A  summary  of  the  compu¬ 
tational  times  for  various  oases  is  shown  in  Table  III-l. 

It  should  be  pointed  out  that  this  CPU  time  includes  the  computation  of 
the  diffracted  mode  amplitudes,  the  physical  optics  approximation,  as  well  as 
the  energy  check.  In  general  most  of  the  CPU  time  is  used  in  filling  up  the 
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Fig.  3-41:  Normalized  Surface  Field  for  a  Sinusoidal  Surface  (d  =  0.  2X, 

a  =  0. 005X  ana  6  =  0)  for  E  polarization. 
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Normalized  Surface  Field  for  a  Sinusoidal  Surface  (d  =  1.9'., 
a  =  0.  25X  and  0  =  0)  for  E  polarization. 
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-43:  Surface  Field  for  a  Triangular  Surface  (d  =  1.75X,  =  20°, 

@2  =  66°  and  B  =  12.2°),  for  E  polarization. 
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j  Number  N  ' 
1  of  Sampling 
“  Points 

Num.ber  M  ! 
of  Terms  in 
the  Infinite 
Series 

Number  Q 
of  Fourier 
terms  used  to 
represent  the 
periodic  surface 

1  CPU  Time  (sec)! 
.(IBM  360/67) 

I 

8 

3 

1 

34.55 

10 

3 

1 

47.  87 

12 

3 

] 

65.55 

10 

5 

1 

96.39 

8 

2 

10 

49.12 

6 

2 

10 

,  34.58 

20 

5 

15 

399. 22 

12 

8 

1 

168.37 

6 

1 

1 

17.13 

4 

3 

i 

14,49 

18 

6 

1 

199.39 

12 

L _ 1 

3 

1 

10 

109.03 

_1 - ! 

TABLE  ni-l:  Computational  Time  for  E  Polarization 


matrix  and  performing  the  matrix  inversion.  For  the  former  operation,  the 

2  3 

CPU  time  varies  as  N  while  for  the  latter  operation  it  is  proportional  to  N  . 

To  see  exactly  how  the  CPU  time  is  distributed  among  the  various  steps  of  the 


computation,  let  us  consider  a  typical  example  where  N  = 

12,  M  =  3  andQ  =  1 

The  distribution  of  the  CPU  time  is 

1 ,  Compulation  of  matrix  elements 

44.  0  sec. 

2.  Matrix  inversion 

8. 3  sec. 

3 .  Physical  optics  approximation 

5. 2  sec. 

4.  Energy  conservation  check  and  diffracted 

8.  0  sec. 

mode  amplitude  calculations 
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Chapter  IV 

SCATTERED  FIELD  FOR  E  POLARIZATION 

In  Chapter  in  we  were  primarily  cim  5e):ned  \ath  the  presentation  of 
numerical  data  for  the  surface  field.  In  g.'^'»3i’al,  however,  the  far  field  is  of 
more  direct  interest,  and  the  solution  of  t  ;  j  integral  equation  for  the  surface 
field  is  then  only  an  intermediate  step.  i’.  this  chapter  we  will  limit  our 
attention  to  the  scattered  field,  the  scai  Tred  energy  and  some  interesting 
physical  phenomena  known  as  anomalies  associated  with  the  diffracted  energy. 
Comparison  is  also  made  with  the  phy'  i^ial  optics  approximation,  and  an  approx¬ 
imate  technique  adequate  for  predicti;'.,  the  field  backscattered  by  a  sinusoidal 
surface  at  oblique  incidence  is  given 

4. 1  Diffracted  Mode  Amplitudes 

When  a  plane  wave  is  inciderJ  on  a  periodic  surface,  the  scattered  field 
can  be  represented  as  an  angular  spectrum  of  plane  waves,  which  spectrum  is 
discrete  by  virtue  of  the  periodic  nature  of  the  boundary  condition  at  the  surface. 
Each  of  the  infinity  of  waves  m>iking  up  the  spectrum  has  associated  with  it  a 
diffraction  angle  which  may  be  real  or  complex  and  is  determined  by  the  grating 
law.  Whereas  the  amplitude  of  i.he  wave  is  a  function  of  the  profile  size  and 
shape,  and  the  directions  of  incidence  and  diffraction,  the  diffraction  angle 
depends  only  on  the  value  of  d/X  and  the  direction  of  incidence. 

A  finite  number  of  diffracted  waves  represent  propagating  modes  and 
these  are  the  Important  ones  far  from  the  boundary.  The  remaining  modes  are 
evanescent  and  though  these  do  not  serve  to  carry  energy  away  from  the  surface, 
they  do  play  a  vital  role  in  affecting  the  amplitudes  of  the  propagating  modes. 

The  number  of  modes  that  propagate  can  be  determined  from  the  expressions 
+  + 

for  given  in  Eq.  (2.30):  is  real,  the  corresponding  mode  propagates 

without  attenuation,  whereas  if  X  ^  is  pure  imaginary,  the  mode  is  evanescent. 

To  find  the  (complex)  amplitudes  of  the  diffracted  waves  we  proceed  as  follows. 
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Let  y  =  fi(x)  be  the  profile  of  the  surface.  Assumii^  the  incident  plane 
wave  to  be 


i  _  -jk<xsin0-ycos0) 
bt  —  e 
z 


(4.1) 


(see  Eq.  (2.20»>  and  by  invoking  the  periodicity  of  the  surface,  the  scattered 
field  can  be  written  as 

00 


E  =  >  A  e  ^  (4.2) 

z  . . *  m 

m=-oo 


for  y  >maxf(x),  where 


n  2m7r  .  ,  ,  . 

B  =  -r-  +ksine  , 
m  d 


and 


X  =\/k^-/3^  • 


The  field  arising  from  the  currents  induced  on  the  surface  is  given  by 


Jo 


E®(x.y)=^  \  P^\/l+{f'(x')y  K^(x')dx' 


(4.3) 


where  is  as  shown  in  Eq.  (2.24).  In  particular,  this  is  valid  for  y  >maxf(x), 
and  hence,  by  combining  Eqs .  (2.22),  (2.24),  (4.2)  and  (4.3)  we  have 


00 
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from  which  we  obtain 


,  j(p  x'+X  y')  I - - - — 

A  =  —^r~  \  e  \/l+ •[f'(x')}'  K  {x')dx'  .  (4.4) 


rn  2dX 


m 


»0 


th 


Notice  that  the  angle  of  diffraction  0^  for  the  m  mode  is  given  by 


tane  =  j3  /X 
m  m  m 


(4.5) 


which,  after  a  little  mathematical  nanipulation,  can  be  reduced  to  the  grating 
law  formula 


2in0  =  ■—  +  sin0  .  (4.6) 

m  d 

Having  determined  the  current  distributions  K  (x)  it  is  therefore  a  trivial 

z 

matter  to  compute  the  amplitudes  v  f  all  the  diffracted  modes,  both  propagating 
and  evanescent,  and  the  amplitudes  of  the  propagating  waves  for  several  dif¬ 
ferent  values  of  d/X,  a/X  and  0  for  a  sinusoidal  surface  are  given  in  Tables 
IV- 1  through  r\^-4.  As  a  basis  for  comparison,  the  physical  optics  approxima¬ 
tions  are  included.  In  general,  the  physical  optics  values,  with  or  without 
shadowing  included,  tend  to  be  too  small,  and  we  further  note  that  if  shadowing 
is  excluded  (or  is  not  present),  the  physical  optics  estimates  are  the  same  for 
both  polarizations  (see  Appendix  A) . 

TABLE  IV- 1 


Amplitudes  of  Diffracted  Waves  for  a  Sinusoidal  Surface 
(d  =  0.2X,  a  =  O.IX). 
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Amplitudes  of  Diffracted  Wa/es  for  a  Sinusoidal  Surfece 
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TABLE  IV-3 

Aijoplitudes  of  Diffracted  Waves  for  a  Sinusoidal  Surface 
(d  =  0.4X,  a  =  0.2X). 
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TABLE  IV-4 

Amplitudes  of  Diffracted  Waves  for  a  Sinusoidal  Surface 
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4.2  Diffracted  Energy  and  Wood  Anomalies 

4.2. 1  Discussion 

As  mentioned  in  Section  4.1,  when  a  plane  wave  is  incident  on  a  periodic 
surface,  the  scattered  field  forms  a  discrete  spectrum  of  plane  waves,  either 
propagating  or  attenuated  away  from  the  surface.  Only  the  propagating  modes 
can  carry  energy  away  from  the  surface  and  are  responsible  for  the  far  field. 
Since  the  surface  is  perfectly  conducting,  the  diffracted  energy  as  a  whole  must 
equal  that  carried  by  the  incident  fieUi. 

The  problem  of  determining  how  the  incident  energy  is  distributed  among 
the  diffracted  waves  of  various  orders  is  a  relitively  old  one,  and  has  been  the 
subject  of  numerous  theoretical  and  experimental  investigations,  many  of  them 
motivated  by  the  desire  for  more  efficient  optical  diffraction.  This  subject  is 
particularly  important  in  the  design  of  diffraction  gratings  for  use  as  mono¬ 
chromators  in  pectroscopy  because  of  their  superior  properties,  such  as  higher 
powers  of  resolution  and  lower  absorption  of  light,  in  comparison  with  prisms. 

If  we  examine  the  spectrum  of  light  resolved  by  optical  diffraction  gratings, 
we  often  find  rapid  variations  in  the  intensity  of  some  diffracted  modes  over  a 
small  range  of  incidence  angles  or  over  a  narrow  frequency  band.  These 
phenomena  are  termed  the  Wood  anomalies.  Unlike  ghosts  (Rowland  or  Lyman), 
anomalies  are  not  caused  by  errors  in  the  spacing  of  the  rulings,  but  are  due 
to  a  type  of  resonant  interaction  among  the  diffracted  modes  scattered  from 
different  periods.  An  anomaly  may  appear  as  a  bright  or  dark  band  in  an  other¬ 
wise  normal  spectrum. 

Historically,  grating  anomalies  were  first  observed  in  reflection  grating 
spectra  by  Wood  in  1902.  Since  these  effects  could  not  be  explained  by  means 
of  ordinary  grating  theory.  Wood  termed  them  "anomalies".  At  that  time,  how¬ 
ever,  he  reported  that  anomalies  could  occur  only  if  the  incident  field  was  H- 
polarized  (S-type  anomalies),  and  there  were  no  anomalies  for  E  polarization 
(P-type  anomalies).  Later,  experiments  by  Ingersoll  (1920)  and  Strong  (1936) 
also  showed  no  evidence  of  P-type  anomalies  in  grating  measurements. 
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The  first  theoretical  treatment  was  given  by  Rayleigh  (1907).  His 
approach  was  based  on  an  expansion  of  the  scattered  field  in  terms  of  outgoing 
waves  only.  (Details  of  this  will  be  given  in  a  later  chapter. )  Using  this 
assumption,  he  found  that  the  scattered  field  was  singular  at  wavelengths  for 
which  one  of  the  diffracted  modes  emerges  from  the  grating  at  the  grazing 
angle.  These  wavelengths  are  now  known  as  the  "Rayleigh  wavelengths", 

^  (  +  l-sin0)  m  =  1,2,3 .  (4.7) 

He  also  observed  that  the  occurrence  of  such  singularities  corresponded  to  the 
appearance  of  the  Wood  anomalies.  Furthermore,  these  singularities  appear 
only  when  the  incident  field  is  H-polarized,  thereby  accounting  for  the  S 
anomalies,  and  if  the  incident  field  is  E  polarized,  the  theory  predicts  a  regular 
behavior  near 

It  should  be  noted  that  although  Rayleigh's  theory  correctly  predicts  the 
major  features  experimentally  observed  at  that  time,  it  does  not  give  the  shupe 
of  the  bands  associated  with  the  S  anomalies  because  of  the  indicated  singularity 
at  the  Rayleigh  wavelength. 

In  his  later  publications.  Wood  (1935)  suggested  the  existence  of  the  P 
anomalies,  but  it  was  Palmer  (1952)  who  succeeded  in  detecting  these  anomalies 
experimentally ,  thus  proving  their  existence  beyond  doubt.  Furthermore, 
Palmer  pointed  out  that  the  P  anomalies  (unlike  the  S-type)  occur  only  for 
gratings  with  deep  grooves.  Since  Rayleigh's  assumption  (chapter  Vn)  is  valid 
for  shallow  grooves  only,  no  inconsistency  is  present,  but  Rayleigh's  theory  is 
incomplete.  Recent,  more  sophisticated  theoretical  treatments  by  Hessel  and 
Oliner  (1965),  Wirgin  (1969),  Itoh  and  Mlttra(1969)  and  many  others  also  indi¬ 
cate  the  existence  of  P  anomalies. 

The  numerical  results  which  we  have  obtained  can  also  be  used  to 
explore  the  Wood  anomalies,  and  this  we  shall  now  do. 
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4.2.2  Computed  Data  for  the  Diffracted  Energy 
a)  Two  Radiating  Modes.  Figure  4-1  shows  the  plot  of  normalized 
diffracted  energy  versus  surface  height  for  a  full-wave  rectified  profile  having 
d  =  0.  75X  with  9  =  41. 8^.  Since  an  equivalent  plot  for  a  sinusoidal  surface 
was  obtained  by  Zaki  (1969),  it  is  included  here  for  comparison  purposes.  For 
this  particular  choice  of  parameters,  the  diffracted  modes  corresponding  to 
m  =  0  and  m  =  -1  can  propagate. 

Examination  of  the  curves  show  that  the  diffracted  energy  for  either  mode 
displays  a  standing  wave  behavior.  At  a  certain  value  of  a/X  depending  on  the 
profile  shape,  there  is  a  complete  exchange  of  energy  between  the  two  modes, 
indicating  the  existence  of  a  P-type  anomaly.  When  this  happens,  the  total 
diffracted  energy  is  carried  by  the  m  =  -1  mode,  and  thus  there  will  be  a  bright 
band  in  the  diffraction  grating  spectrum  for  the  m  =  -1  mode.  The  value  of  a/X 
for  which  it  occurs  is  (Fig.  4-1)  0.  8  for  a  sinusoidal  surface,  but  just  over  2.4 
for  the  full-wave  rectified  surface.  Note,  however,  that  the  actual  peak-to- 
trough  depth  for  the  sinusoidal  surface  is  2a,  so  that  as  a  function  of  the  total 
depth  the  change  is  not  quite  so  dramatic. 

The  above  data  clearly  demonstrate  the  P  anomalies  are  deep-grooved 
phenomena  and  dependent  on  the  shape  of  the  grooves.  These  are  consistent 
with  Palmer's  experimental  observations.  Furthermore,  the  curves  show  that 
the  band  for  the  full-wave  rectified  profile  is  broader  than  that  for  the  sinusoidal 
surface,  and  hence  for  the  former  surface  the  diffracted  energy  in  both  modes 
is  less  sensitive  to  a  change  in  surface  height. 

Note  that  these  P  anomalies  differ  from  the  classical  S  anomalies  in 
that  they  are  broader  and  need  not  occur  near  the  Rayleigh  wavelengths.  Under 
certain  conditions,  a  total  (ICO  per  cent)  enei*gy  conversion  takes  place  between 
the  incident  plane  wave  and  tlie  m  =  -1  mode,  with  the  amplitude  of  the  specularly- 
reflected  wave  becoming  zero  at  this  point.  These  phenomena  were  also  investi¬ 
gated  by  Tseng,  Hessei  and  Oliner  (1968),  and  tliey  termed  .Jhem  Brewster-angle 
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effects  in  an  analogue  to  the  well-known  Brewster-angle  effect  for  the  case  of  an 
interface  between  two  different,  homogeneous,  isotropic  media.  Furthermore, 
they  pointed  out  that  at  the  first  P  anomaly  we  must  have  the  following  condition 
(Bragg  condition): 


kdsinO  •- TT  (4.8) 

Therefore,  if  we  substitute  this  relation  into  Eq.  (2,21),  we  obtain 

K  (r+rndx)  =  K  (r)e  m  =  tl,  1 2, . . . 

z  —  2  — 

("K  (r<  if  m  =  even  (4,9) 

I  -K  (r)  =  odd 

^  z  — 


which  shows  that  at  this  anomaly  the  current  distribution  over  any  chosen  cell 
(or  period)  must  be  in  or  out  of  phase  with  that  over  the  reference  cell  (m  =  0) 
depending  on  whether  the  chosen  cell  is  displaced  an  even  or  odd  number  of 
periods,  respectively,  from  the  reference  cell.  But  this  is  only  a  necessary 
condition  for  the  existence  of  an  anomaly.  There  does  not  appear  to  be  any 
rule  for  predicting  the  exact  depth  of  groove  for  an  anomaly  to  occur  since  this 
depth  is  different  for  different  profiles.  However,  using  the  condition  (4.8),  we 
can  derive  some  information  about  the  nature  of  the  current  distribution  when 
the  first  anomaly  occurs.  This  is  done  as  follows. 

From  Eqs.  (4.2)  and  (4.8),  for  the  m  =  0  mode: 


=  k  sin  9  =  TT /d  , 
=  kcosQ  ; 


(4. 10) 


and  for  the  m  =  -1  mode: 

=  -Tf/d  , 

X  j  =  -kcos9  . 


(4.11) 
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Substituting  Eq.  (4.10)  into  Eq.  (4.4)  we  get 


C" 

A  =  -Ji’iL—  \ 

0  idkcosO  \  ® 

Jo 


^  j(jx+kcos0y), - ^ 

e  \/l  +  f'(x)  K  (x)dx. 

z 


(4. 12) 


Similarly,  using  Eq.  (4.11), 


\d  -j(--x+kcos0y) 


A  —  1 

-l~2dkcos6  1  ® 


/l  +  f'(x)^  K  (x)dx  ,  (4.13) 


which  can  be  written  as 


,  -f  X  j(^x-^kcos8yV.- 

-l  2dkcos0  I  ®  c  J 

JO 


'r - 2 

JH-f'(x)  K  (x)dx  .  (4.14) 

■  r? 


Furthermore,  at  the  first  P  anomaly: 


'W  =  “  • 


(4.15) 


Observe  that  Eqs.  (4. 12)  and  (4. 14)  differ  only  by  a  factor  exp(-j  -~x) 
m  the  kernel,  but  their  magnitudes  differ  by  unity.  To  show  that  this  is  also 
true  if  we  replace  exp(-j  ~  x)  by  its  conjugate  (exp(j  ^  x^  ,  change  the 
variable  x  to  -x  to  get 


j("Tx+kyco8  0) 

— a  \  S  \'l+f'(-x)  K  (-x)dx. 

2dkcos0  S  z 

Jo 


(4.16) 


If  it  is  assumed  that  f(x)  is  an  even  function  then 


j(-^  x  +  kyco80) 

A  H  \  ®  '^l+f'(x)  K  (-x)dx  , 

- 1  2dk  cos  0  1  *  z 


(4.17) 


9'^, 


and  using  Eq.  (2.21) 


=  K  (x)  e 
z 


.  27r 


(4.18) 


Thus  Eq,  (4.17)  becomes 


^1  = 


(j/i 


2dk  cos  6 


\d  .  27r  ,/  71^  .  ,  a  \ 

J-J-x  j(-x+kcosey) 

e  e 


(/l  +  f'(x)' 


K  (x)dx 
z 


(4.19) 


which  differs  from  Eq.  (4. 14)  only  in  sign.  Moreover,  if  we  add  Eqs.  (4, 14) 
and  (4. 19),  we  have 


dk  cos  6 


j('^  x+kcosey)  j - 2 

!  \i'l+f'(x)  K  (x)dx  , 

•  z 


(4.20) 


whereas  if  we  subtract  Eq.  (4. 19)  from  Eq.  (4. 14), 

id 


0  = 


_wiL 


dk  cos  6 


i'2ir  '  K-^Jx+kcosey)  r . 


cos  t  —  x  e 


'Vd 


Vl  +  f'(x)  K  (x)dx  (4.21) 
z 


lo 


Therefore  from  Eqs.  (4.12),  (4,20)  and  (4.21)  and  the  orthogonality 
property  of  the  trigonometric  function,  we  can  conclude  that 


K  (x)  = 
z 


kd  COB  9 


\’l  +  f'(x)^ 


-j(  x+kcos  0y) 
e 


(4.22) 


which  describes  the  behavior  of  the  surface  field  at  the  first  P  anomaly.  In 
particular,  if  we  consider  a  sinusoidal  surface,  then  the  modulus  of  the  surface 
field  will  exhibit  a  standing  wave  pattern  showing  the  propagation  of  unattenuated 
surface  waves  along  the  x-axls.  Once  again  this  phenomenon  is  seen  to  be 
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analogous  to  the  Brewster  angle  effect  which  is  used  in  launching  a  surface 
v/ave  along  an  interfece  of  two  homogeneous,  isotropic  media.  For  a  perfectly 
conducting  flat  sheet,  this  Brewster  angle  is  of  course  equal  *o  90°  (with  the 
normal  to  the  surface),  but  since  the  surface  is  corrugated,  it  can  be  modelled 
by  an  appropriate  reactance  surface  (see,  for  exam.ple,  Hessei  and  Oliner,  1965), 
whose  Brewster  angle  can  now  assume  values  other  than  90°, 

Another  illustration  for  two  radiating  modes  is  shown  in  Fig.  4-2  for  a 
full-wave  rectified  profile  at  oblique  incidence  with  9  =  60°.  The  parameters 
of  the  profile  are  a  *  0, d  =  0. 30\  and  the  period  varies  from  0  to  0. 7X. 

The  plots  show  that  there  is  a  gradual  exchange  of  diffracted  energy  between  the 
m  =  0  and  m  =  -1  modes  when  d  >  0.535X,  which  is  a  Rayleigh  wavelength. 

b)  Three  Radiating  Modes.  As  the  period  becomes  larger,  more  radiating 
modes  are  possible.  A  case  of  three  radiating  modes  is  shown  in  Fig.  4-3  for  a 
sinusoidal  surface  at  oblique  Incidence,  6  =  60°.  The  parameters  of  the  surface 
are  d  =  1.1 55X,  and  a  varies  from  0.05  to  0.7X.  The  three  radiating  modes 
are  the  m  =  0,  m  =  -1  and  m  =  -2  modes,  and  as  computed  from  the  grating 
law,  their  diffracted  angles  are  60°  (m  =  0),  0  (m  =  -1)  and  -60°  (m  =  -2). 
Thus,  the  m  =  0  mode  is  the  specular  mode,  radiating  in  the  foreward  scattered 
direction,  the  m  =  -1  mode  radiates  in  the  normal  direction,  and  the  m  =  -2 
mode  is  the  backscattered  mode,  radiating  in  the  reverse  direction  of  incidence. 

The  plots  show  a  gradual  exchange  of  the  diffracted  energy  among  the 
three  radiating  modes.  For  a  <  0.05X,  the  dominant  mode  is  the  m  =  0  mode, 
which  carries  almost  all  the  diffracted  energy.  Thus  from  the  far  field  point 
of  view,  the  surface  is  similar  to  a  flat  surface.  As  the  amplitude  increases, 
the  other  two  modes  come  in,  and  play  an  important  role  in  sharing  the  radiated 
energy.  Since  no  single  mode  is  observed  to  carry  the  whole  (or  none)  of  the 
energy  (except,  of  course,  when  a«0.05X),  even  for  a  as  large  as  ().7X, 
it  is  unlikely  that  the  P  a..omalies  will  occur  with  three  radiating  modes  when 
the  surface  height  is  not  very  deep  (of  order  X  or  less);  and  if  they  do  occur,  it 
will  only  be  for  groove  depths  much  greater  than  those  for  two  radiating  modes. 
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Fig.  4-3:  Diffracted  Energy  for  a  Sinusoidal  Surface  (d  =  1, 155X  and 

6  =  60°)  for  E  polarization. 
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The  same  kind  of  plots  were  also  obtained  using  the  physical  optics 
approximation,  either  with  or  without  shadowing,  and  they  are  shown  in  Figs. 
4-4  and  4-5.  Comparison  with  Fig.  1-3  indicates  that  neither  version  of  the 
physical  optics  method  gives  satisfactory  agreement  with  the  exact  solution 
except  for  the  range  a  <  0.  IX.  For  such  surface  heights,  there  is  no  shadow 
region. 

4.3  Conservation  of  Energy 

One  criterion  that  we  have  used  to  check  the  accuracy  of  the  numerical 
solution  is  conservation  of  energy.  The  numerical  results  do  indeed  satisfy 
this  condition,  as  can  be  seen  from  Figs.  4-1  and  4-2  for  the  two  radiating  mode 
case,  and  from  Fig.  4-3  for  the  three  radiating  mode  case.  The  results  for  a 
single  radiating  mode  are  given  in  Table  IV-5. 


TABLE  IV-5 

Conservation  of  Energy  for  a  Sinusoidal  Surface 
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0 
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60 

1 
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0.25 

0 
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i 
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-57.98 
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. J 

From  Table  IV-5  and  Figs.  4-4  and  4-5,  we  see  that  the  physical  optics 
approximations  In  general  do  not  satisfy  the  condition  of  energy  conservation. 

Another  plot  showing  the  relation  between  the  diffracted  energy  and  the 
minimum  radius  of  curvature  (at  the  surface  peak  or  trough)  is  shown  in 
Fig.  4-6.  It  indicates  that  for  >4.5,  most  of  the  diffracted  energy  is 
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Surface  Amplitude  (a/X) 


Fig.  4-5:  Diffracted  Energy  for  a  Sinusoidal  Surface  (d  =  1. 155X  and 

0  =  60  ),  Computed  Using  the  Modified  Physical  Optics 
Method  (without  shadow). 
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C£2i"i0d  away  by  the  m  =  0  mode  so  that  as  far  as  the  radiated  field  is  concerned, 

the  sinusoidal  surface  can  be  considered  a  flat  surface  if  kp  ,  >4.5. 

min 

4.4  Approximate  Techniques 

Although  the  main  purpose  of  this  work  is  to  develop  numerical  procedures 
for  solving  the  integral  equations  for  scattering  Jh:om  periodic  surfaces,  a  sub¬ 
sidiary  objective  is  to  use  the  understsinding  gained  from  our  results  to  treat  back- 
scattering  from  rough  surfaces,  with  particular  reference  to  oblique  angles  of 
incidence.  A  desirable  step  in  this  process  is  to  find  more  simple  methods  for 
predicting  the  field  back-scattered  by  a  periodic  surface  without  resorting  to 
the  integral  equation. 

The  method  that  naturally  comes  to  mind  is  the  physical  optics  method,  but 
as  we  have  already  seen,  the  predictions  are  rather  poor  in  the  far  field  and  even 
worse  as  regards  the  near  field.  One  reason  for  this  is  that  the  aq)proxlmation 
is  a  local  one  onjly,  taking  no  account  of  the  coupling  between  cells  and  of  the 
interactions  among  the  diffracted  energy. 

If  we  examine  the  curves  shown  in  Figs.  3-26  through  3-30  for  the 
sinusoidal  surface,  we  notice  that  over  the  illuminated  portion  of  the  surface 
and  particularly  near  the  surface  peak,  the  current  is  primarily  determined  by 
the  local  profile  shape,  and  is  affected  relatively  little  by  the  other  part  of  the 
surface.  Furthermore,  in  most  cases  the  surface  field  is  largest  over  this 
region,  and  the  phase  does  not  show  the  rapid  variation  which  it  does  in  the 
shadow  region.  These  facts  reinforced  with  some  reasoning  based  on  ray  optics 
suggest  that  it  may  be  possible  to  estimate  the  back  scattered  field  with  reason¬ 
able  accuracy  using  only  a  knowledge  of  the  current  in  the  illuminated  region. 

This  is  indeed  the  case,  as  was  shown  by  applying  the  stationary  phase  method 
to  Eq,  (4.4)  over  the  illuminated  region  and  using  the  exact  values  for  the  surface 
field.  The  results  are  presented  in  Table  lV-6  and  compared  with  the  exact 
results. 
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TABLE  IV-  6 


r 

I 


I 

t 


L 


Comparison  of  the  energy  backscatterod.  by  a  sinusoidal  surface: 
d  =  1 . 155X,  0  =  69°,  computed  using  the  exact  and  stationary 
phase  methods. 


Surface  i 

Amplitude  (a/X)  i 


Normalized  Backscattered  Energy 
Exact  Method  Stationary  Phase 


0.3 

0.176 

0.152 

0.4 

0.280 

0.320 

0.5 

0.320 

0.360 

0.6 

0.272 

0.344 

0.7 

1 

0.200 

0.216 

Since  the  above  computation  requires  a  knowledge  of  the  true  surface 
field,  albeit  over  only  a  limited  portion  of  the  surface,  we  cannot  avoid  solving 
the  integral  equation  unless  we  can  find  an  approximate  method  for  estimating 
the  surface  field  in  this  region.  Because  of  the  relatively  good  agreement  be¬ 
tween  the  physical  optics  phase  and  the  exact  phase  over  the  illuminateo  region, 
the  phase  in  this  region  can  be  simply  taken  as  the  incident  phase.  As  for  the 
current  modulus,  we  can  approximate  it  in  either  of  two  different  ways . 
a.  Circular  cylinder  approximation. 

For  an  E-polarized  plane  wave  liicident  on  a  circular  cylinder  of  radius 
a,  the  modulus  of  the  current  induced  at.  the  point  P(x,  y)  is 


mod|k  (x)j  = 


00 


i4iira 


n=-oo 


.-a  jnu(x) 
_ 

i/'''^{ka) 


(4.23) 


where  u(x)  =  obtuse  angle  between  the  viormal  at  P(x,y)  and  the  incident  ray. 
Let  us  now  identify  a  with  the  local  rad'us  of  curvature,  p(x),  at  the  point 
P(x, y)  on  the  surface  (see  Fig.  4-7).  For  a  cdnufioidal  surface,  p(x)  is  given 
by 


I 


Fig.  4-7; 


Approximation  of  the  Surface  Current  Modulus  Over  the 
Illuminated  Region  of  the  Sinusoidal  Surface  by  that  on  a 
Circular  Cylinder. 
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This  equation  can  now  be  used  to  compute  ^he  current  modulus  in  the  illuminated 
region  of  the  surface. 

b.  Parabolic  cylinder  approximation. 

Another  way  of  approximating  the  current  modulus  over  the  lit  region  of 
the  sinusoidal  surface  is  to  use  our  knowledge  of  the  current  induced  on  a  per¬ 
fectly  conducting  parabolic  ^iyiinder  whose  vertex  coincides  with  the  peak  of 
the  sinusoidal  surface  (se*?;  Fig.  4-8),  and  both  surfaces  have  the  same  radius 
of  curvature  at  this  point.  Because  of  the  mathematical  complexity  of  the  ex¬ 
pression  when  the  cylinder  is  at  oblique  incidence,  we  use  instead  the  normal 
incidence  results  and  modify  it  by  the  factor  cosu  defined  above. 

At  nonnal.  incidence ,  the  current  induced  on  a  parabolic  cylinder  is 


K  (x)  - 
z 


— 


F(n^) 


1  „  2, 
--jk(?  +r)  ) 


(4.27) 


where  ?,  r),  z  are  parabolic  cylindrical  coordinates  related  to  the  rectangular 
Cartesian  coordinates  by  the  transformation 


Approximation  of  the  Surface  Current  Modulus 
over  the  Illuminated  Region  of  the  Sinusoidal  Surface 
by  that  on  a  Parabolic  Cylinder, 
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-  i  /c2  2,  _ 

x=-(?-n),  y  =  T)?,  2  =  z  , 


(4.28) 


and  F(TjVh)  is  the  Fresnel  integral  which  can  be  computed  using 


Thus,  multiplying  Eq.  (4.27)  by  cosu,  we  get 


(4.29) 


1 

F(r)\flc) 


1  •i./(»2  1 

-gjW?  +J7  ) 
e  cosu 


(4.31) 


which  is  the  equation  used  to  compute  the  current  modulus  in  the  illuminated 
region  of  a  sinusoidal  surface  at  oblique  incidence. 

Frcir  experience,  it  has  been  found  that  the  parabolic  cylinder 
approximate  ■' ;  generally  gives  results  slightly  less  than  the  exact  values, 
whereas  the  circular  cylinder  approximation  gives  slightly  larger  values. 

The  rebults  are  shown  in  Table  IV-7,  and  also  compared  with  the  exact  values. 


TABLE  IV-7 

Comparison  of  Normalized  Current  Modulus  over  the  Illuminated 
Region  for  the  Sinusoidal  Surfece  with  d  =  1.155X  and  $  =  60°. 

I  Location  of  Point  T  rjfrcular  ]'  Paribblic) 

Amplitude!  On  the  Surface,  Exact  Cylinder  Cylinder  i 
_ ^  Results  j  Approx.  '>  Approx.  | 


Chapter  V 

SURFACE  FIELD  DATA  FOR  H  POLARIZATION 

In  this  chapter  surface  field  distributions  obtained  by  numerical  solution 
of  the  integral  equation  for  H  polarization  are  presented.  The  particular 
profiles  considered  are  sinusoidal,  full-wave  rectified,  inverted  full-wave 
rectified,  and  triangular  (see  Fig.  3-1)  with  various  periods  and  angles  of 
incidence.  In  each  case  the  physical  optics  approximation  is  presented  as  a 
basis  for  comparison,  and  comparison  also  made  with  analogous  results 
for  E  polarization.  Many  of  the  remax’ks  made  in  Chapter  ni  for  E  polarization 
hold  for  H  polarization  as  well,  but  the  surface  field  now  is  generally  more 
complicated  and  unpredictable  even  over  a  small  region  of  the  surface. 

According  to  physical  optics,  the  surface  current  induced  on  the  surface 
is 

(5.1) 

in  the  Illuminated  region,  and 
K^(x)  =  0 

in  the  shadow. 

As  in  the  case  of  E  polarization,  however,  this  approximation  is  inadequate 
for  many  purposes. 

5.1  Surfaces  of  Small  Period  (d/.X^  0.2) 

Figures  5-1  through  5-12  represent  the  surface  fields  on  periodic  surfaces 
of  four  different  profiles,  each  having  a  period  d  =  0.2\  and  a  maximum  ampli¬ 
tude  a  =  0.  IX,  The  angles  of  incidence  are  0  =  0,  30*^,  and  60°. 

Let  us  first  examine  the  results  for  a  sinusoidal  surface  (Figs.  5-1  through 
5-3).  The  current  modulus  now  has  a  maximum  at  the  surface  trough  and  a 
minimum  at  the  peak,  and  in  contrast  to  the  case  for  E  polarization,  it  remains 
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Phase  (degrees)  Normalized  Current  Modulus 


Ill 


substantially  unchanged  as  the  angle  of  incidence  is  varied.  The  physical  optics 
approximation  is  somewhat  like  the  mean  of  the  exact  solution.  For  0  =  0,  the 
phase  is  almost  constant  (as  for  E  polarization),  but  there  is  more  variation, 
especially  near  the  surface  peak,  when  0^0. 

Figures  5-4  through  5-9  show  the  surface  fields  for  full-wave  rectified 
and  inverted  full-wave  rectified  profiles.  Except  that  there  is  a  slight  improve¬ 
ment  in  the  agreement  between  the  exact  modulus  and  the  modulus  of  the  physical 
optics,  the  general  behavior  is  rather  similar  to  that  for  a  sinusoidal  surface. 
Shadowing  has  little  effect  on  the  current  modulus. 

The  last  three  figures  (Figs.  5-10  through  5-12)  are  for  triangular 
profiles,  and  we  note  only  the  absence  of  any  "edge  effects". 

A  point  worth  mentioning  is  that  for  all  the  profiles  caisidered,  the 
exact  phase  at  normal  incidence  is  close  (in  a  mean  sense  at  least)  to  the 
physical  optics  phase.  In  order  to  see  why  the  modulus  of  the  surface  field  is  so 
insensitive  to  any  change  in  the  angle  of  incidence,  let  us  go  back  to  Eq.  (2,35). 
As  in  the  case  of  E  polarization,  if  the  period  is  sufficiently  small,  we  can 
approximate  the  modified  Green's  function,  Eq.  (2.36),  by  the  first  two  and 
the  last  two  terms.  Thus,  the  approximate  integral  equation  becomes 

5G2(x,y;x',y')K^(x')dx'2:  jd{2H^(x)-K^(x)]  ,  (5.2) 

0 

where 

G2(x,yjx',y')  = -j  sgn(y-y')-j  |^f'(x')tan0-sgn(y-y')| 

-jksin0(x-x')  -  jkco80[y-y'| 

0  • 

+  ^  ■^f'(x')+  jsgn(y-y')|  cotj^-  ^(x-x')+j  |y-y'j|  + 

■[-f'(x')+j8gn(y-y')j- cot-|^~(x-x')+j b'y'!}"  • 


Normalized  Surface  Field  for  an  Inverted  Full-Wave  Rectified 
Surface  (d  =  0. 2X,  a  =  0.  IX  and  6  -  60°)  for  H  polarization, 

— 0—  exact,  —  physical  optics. 


— ».]  Shadow 


—  Shadow 


Fig.  5-9: 


Normalized  Surface  Field  for  a  Full-Wave  Rectified  Surface 
(d  =  0. 2X,  a  =  0.  IX  and  Q  =  60°)  for  H  polarization,  — o — 
exact,  —  physical  optics. 


Phase  (de^ees)  Normalized  Current  Modulus 
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Fig.  5-10:  Normalized  Surface  Field  for  a  Triar^alar  Surface  (d  =  0.2X., 
f  =  0.  IX,  a  =  0.  IX  and  0  =  0)  for  H  polarization,  — o —  exact, 
—  physical  optics. 
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Fig.  5-11;  Normalized  Surface  Field  for  a  Triangular  Surface 

(d  =  0. 2A.,  i  =  0.  IX,  a  =  0.  IX  and  6  =  30^)  for  H  polarization, 
—0 —  exact,  —  physical  optics. 


Phase  (degrees)  ^  |  Normalised  Current  Modulus 


Fig,  5-12:  Normalized  Surface  Field  for  a  Triangular  Surface  (d  =  0.2X, 
i  -  0. 1\,  a  =  0,  IX  and  6  =  60  )  for  H  polarization  — o~  exact, 
—  physical  optics. 


121 


Note  that  of  the  fotir  terms  in  G^,  only  the  second  contains  the  angle 
of  incidence  explicitly.  Since  the  integral  ^nation  is  singular,  a  rough  picture 
of  the  solution  may  be  obtained  if  we  consider  the  integral  equation  in  the  neigh¬ 
borhood  of  the  singularity.  Thus 

lim  G„(x,  y;  x’,  y*)  =  -j  sgniy  -  y')  -  j  ff'(x*)tane-  sgn(y  -  y’)'l-  + 


x-x 

jy-y' 


— ^0 

0 


•^{f'(x)  +  jsgn(y-y')\  lim  cot  f- <x-x')+j  4  |y-y’ll  + 

^  ^!x-x*!-^0  ^  ^  a  'J 

ly-y’Ko 

■2^f(x')'*'Jsgn(y-y')}  linj  cot  ft  (x-x')+j  ly-y']!.  (5.4) 

lx-x'|~>0  ^  J 

jy-  y']— >0 


Using  the  relation 


cot(s)  =  j 


we  obtain  for  the  cotangent  terms 

2b+j2a 


oot(-a+jb)  =  j 


1  +e 


1  -  e 


2b+j2a 


,  and  cot{a+jb)  =  j 


1+e 


2b-j2a 


1  -  e 


2b-j2a  ' 


(5.5) 


where  a  =  ^(x-x')  ,  b  =  -jly-y')  . 


Therefore, 

lim  cot-j  -  •^(x-x')+j  y-y'll  =  lim 

|x-x'|-~»0  ^  °  °  J  |ai-^0  1- 

ly-y'h»o 


+  e 


2b+j2a) 


|a|->C 
b— >0 


2b+j2a 


2,  liin  z(£L±M 

Axiu  9  9  * 

)al~>0  a  +b 
b-^>  0 


(5.6) 
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Similarly,  for  the  other  term 


lim 

Ix-x*  |-»0 

jy-y']—^ 


(x-x*)+j 


2;  lim 
|a|— >0  a  +b 
b“*>0 


(5,7) 


Substituting  Eqs.  (5.6)  and  (5.7)  into  Eq,  (5.4),  we  obtain,  after  some  mani¬ 
pulation. 


lim  G  (x,y;x',y‘)  = -jsgn(y-y')-jrr(x*)tan0-sgn(y-y')l  + 
lx-x«l->0  ^  1-  J 

jy-y'l— >0 

+  ^ — 2  [b8gn(y-y')+af'(x'^ 

|ap>0  d(a  +b  ) 
b  — >0 

Thus,  in  the  neighborhood  of  the  singularity,  Eq.  (5.2)  becomes 

pX  +  € 

\  •^-jsgn(y~y')- j  jf'(x')tan0-sgn(y-y*^  + 

Jx-e 


(5.8) 


+  ,  lim  — ;; — -  [b8gn(y-y«)->-af'(ji]\  K  (x')dx' 

|ah»0  d(a‘'+b^)  i 

b  — >0 

^  Jd{2H^(x)-K^(x)]  (5.9) 


where  [x-x'l  <e  with  e  (>0)  an  arbitrarily  small  number. 

If  we  assume  K^(x)  to  remain  constant  over  this  small  region  of 
integration,  then  by  the  property  of  sgn(y-y'>,  all  integrations  involving 
3gn(y  -  y')  vanish.  Hence 


AX+e 


2Jaf'(x') 


^-jf'(x')tan0+  lim 

|a|“->0  d(a  +b  )  J 
b  — >0 


K^(x')dx'  rir  jd{2H^^(x)  -  K^(x)  j 


(5.10) 
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If  the  surface  amplitude  is  relatively  small,  a  and  b  can  be  assumed 
to  be  of  the  same  order  of  magnitude,  in  which  case  Eq.  (5.10)  becomes 


D-jf  (x')tan0  +  -  \  K>-)dx'  2^  jd(2H^.(x)  -  K.(x)l  . 

7r(x-x‘)i  ‘  (.  z  t  3 


(5.11) 


which  can  be  solved  by  treating  K  (x)  and  f(x)  as  constant  over  the  range  of 
integration.  We  thus  obtain 


-jf(x)tan0  K^(x)2e  jd  ■j^2H^(x) -  K^(x)|  , 


(5.12) 


and  since  e  is  an  arbitrarily  small  number,  the  left  hand  side  can  be  neglected 
in  comparison  with  the  second  term  on  the  right  hand  side.  Equation  (5.12)  now 
reduces  to 


K.(x)=2H^x),  (5.13) 

V  Z 

which  is  just  the  geometrical  optics  value.  This  is  consistent  with  the  numerical 
results. 

For  a  sinusoidal  profile,  the  effect  on  the  surface  field  produced  by  a 
change  in  the  surface  amplitude  as  the  period  is  kept  constant  is  shown  in  Figs. 
5-1,  5-13  and  5-14.  As  the  surface  height  is  increased,  the  numerical  results 
approach  those  for  a  flat  sheet. 

5.2  Surfaces  of  Slightly  Larger  Period  (0.2.^d/\^  0. 6) 

Figures  5-15  and  5-16  show  the  surface  field  for  a  sinusoidal  profile 
with  d  =  0.4X,  a  =  0.2\  and  0  =  0,  60*^  respectively.  By  comparing  these 
with  Figs.  5-1  and  5-3  it  is  seen  that  doubling  both  d  and  a  has  a  somewhat 
larger  effect  on  the  current  modulus  than  was  the  case  for  E  polarization.  In 
particular,  there  is  a  noticeable  change  of  shape,  with  minimum  no  longer 
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occurring  at  the  surface  peak.  Ouce  again,  there  is  no  agreement  between 
the  exact  data  and  physical  optics,  either  in  modulus  or  phase. 

Results  for  a  full-wave  rectified  profile  with  d  =  0.6\,  a  =  0.3A.  and 
0  =  60*^  are  given  in  Fig.  5-17.  The  undulations  of  the  modulus  clearly  indi¬ 
cate  strong  coupling  among  cells  (or  multiple  scattering  among  ce”~).  Most 
of  the  surface  current  concentra*ec  in  the  illuminated  region.  As  the  angle  of 
incidence  is  further  increased  to  0  *  75^  (Fig.  5-18),  the  modulus  undergoes 
a  marked  change,  whereas  the  phase  remains  constant  over  most  of  the  surface. 

5.3  Surfaces  of  Larger  Period  (d/A  «b0.6) 

Surface  field  distributions  for  surfaces  of  large  period  are  generally 
more  complicated  than  those  for  surfaces  of  small  period,  primarily  because 
of  the  greater  number  of  radiating  modes . 

Let  us  first  consider  the  results  for  two  sinusoidal  surfaces  whose 
periods  are  near  the  Rayleigh  wavelength.  The  first  one  (Fig.  5-19)  has 
d  =  0. 95\,  a  =  0. 25X  and  0=0.  It  is  seen  that  the  current  modulus  curve 
has  a  maximum  at  the  surface  trough  and  minima  at  the  places  where  the 
surface  crosses  the  x  axis.  The  phase,  however,  remains  fairly  constant  over 
the  lower  part  of  the  concave  portion  of  the  surface,  but  changes  rapidly  from 
positive  to  negative  where  the  modulus  has  its  minimum.  The  second  one  is 
given  in  Fig.  5-20,  where  d  =  1.90X,  a  =  0.25X  and  0  =  0.  Here  the  modulus 
has  a  marked  oscillation,  quite  distinct  from  that  found  for  E  polarization.  One 
of  the  maxima  is  at  the  surface  trough,  and  the  current  is  again  a  minimum  at 
the  surface  peak.  The  agreement  with  physical  optics  is  poorer  than  for  E 
polarization,  with  neither  modulus  nor  phase  being  approximated  to  any  real 
extent. 

A  careful  study  of  Figs.  5-19,  5-20  and  the  corresponding  ones  for  E 
polarization  (Figs.  3-20  and  3-21)  leads  us  to  conjecture  that  there  are  strong 
surface  waves  propagating  along  the  x-axis  when  d  is  close  to  2X  for  H 
polarization,  and  when  d  is  close  to  the  wavelength  for  E  polarization. 


Fig.  5-17:  Normalized  Surface  Field  for  a  Full-Wave  Rectified  Surface 

(d  =  0. 6X,  a  =  0. 3X  and  9  =  60°)  for  H  polarization,  — o —  exact, 
—  physical  optics. 
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—  physical  optics. 


Fig.  5-20:  Normalized  Surface  Field  for  a  Sinusoidal  Surface 
(d  =  1.  9X,  a  =  0. 25X  and  0  =  0)  for  H  polarization. 
— o—  exact,  —  physical  optics. 
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Another  case  also  showing  a  marked  oscillation  in  the  modulus  is  given 

in  Fig.  5-21  where  the  surface  is  a  sinusoidal  with  d  =  1.6X,  a  =  0,25X  and 

0  =  45°  (kp  -  =1.63).  The  physical  optics  phase  is  in  better  agreement  wiCh 
min 

the  exact  lAase  except  in  the  region  where  ihe  surface  is  most  likely  to  receive 
shadowing. 

T  o  examine  the  effect  on  the  surface  field  due  to  a  gradual  increase  in 
the  surface  height  when  the  period  and  angle  of  incidence  are  kept  fixed^  we 
select  a  sinusoidal  profile  with  d  =  1. 155\,  0  =  60°  and  allow  a  to  increase 
from  0.05\  to  0.7X.  The  numerical  results  are  shown  in  ^igs.  5-22  through 
5-29.  Again  as  in  the  case  of  E  polarization,  the  agreement  between  the  exact 
results  and  physical  optics  deteriorates  as  a  is  increased.  Unlike  E  polarization, 
however,  this  agreement  remains  exceptionally  good  for  the  phase  even  when 
the  amplitude  becomes  as  large  as  0.2X.  Furthermore,  the  exact  phase  in  the 
illuminated  region  is  very  close  to  that  of  the  incident  field  in  all  cases  con¬ 
sidered,  but  physical  optics  overestimates  the  modulus. 

To  see  how  much  of  the  surface  field  depends  on  the  concavity  of  the 
surface,  we  replaced  the  concave  portion  (d/4  <  x  <  3d/4)  of  the  sinusoidal 
surface  (a  =  0.5X)  by  a  flat  part.  The  results  are  presented  in  Fig.  5-30. 
Comparison  with  Fig.  5-27  shows  that  this  replacement  has  a  marked  effect  on 
the  surface  field  in  the  shadow  region,  but  little  effect  elsewhere. 

Additional  data  for  a  full-wave  rectified  profile  (d  =  0. 85X  a  =  0. 30X  and 
0  =  36°),  and  a  triangular  profile  (d  =  1.75X,  0^  =  20°,  =  66°,  and 

0  =  12.2°)  are  given  in  Figs.  5-31  and  5-32.  fii  the  former  case,  the  exact 
modulus  and  phase  are  in  good  agreement  with  the  physical  optics  values  over 
the  left  hand  portion  of  the  surface,  but  the  agreement  (particularly  for  the 
modulus)  is  less  good  on  the  right.  For  the  triangular  profile,  the  exact 
modulus  shows  an  oscillatory  variation  not  predicted  by  physical  optics,  and 
this  oscillation  is  most  pronounced  on  the  far  right  hand  portion  of  the  surface. 
The  physical  optics  phase  is  again  a  good  approxinmtion  to  the  exact  phase. 
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Fig  5-29:  Normalized  Surface  Field  for  a  Sinusoidal  Surface  (d  1. 153X, 

a  =  0.7X  and  6  •  60®),  for  H  polarization,  — o--  exact, 

—  physical  optics. 
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Figures  5-33  through  5-37  illustrate  the  effect  on  the  surface  current 
distribution  produced  by  a  progressive  change  in  the  period  of  the  surfece  while 
the  other  parameters  are  kept  fixed.  The  surface  selected  is  a  full-wave 
rectified  surface  with  a  *  0. 3\,  d  =  0. 3A,  0  =  60^  and  f  {the  distance  between 
two  neighboring  periods) varying  from  0. 3\  to  0. 7\  in  0. 1\  steps.  For  the 
first  two  figures,  where  f  ^  0.4\  the  modulus  shows  very  good  symmetry  in 
spite  of  shadowing,  and  over  the  central  portion  of  the  surface  at  least  the 
phase  is  almost  constant.  The  physical  optics  estimates  are  rather  unsatis¬ 
factory.  As  d  becomes  larger  and  approaches  the  Rayleigh  wavelength  (equal 
to  0. 535X,  according  to  Eq.  (4. 7)),  the  surface  field  experiences  a  marked 
change  as  shown  in  Fig.  5-35,  but  the  modulus  still  preserves  its  symmetry. 
Fuither  increase  in  d  again  brings  noticeable  changes  in  the  shape  of  the 
modulus  (see  Figs.  5-36  and  5-37),  and  in  spite  of  the  factldwit  the  radius  of 
curvature  is  a  minimum  at  the  surface  peak,  the  physical  optics  estimate  of 
the  modulus  is  most  accurate  in  that  region.  Compared  v/ith  the  analogous 
results  for  E-polarization,  the  surface  fields  for  these  relatively  small  separ¬ 
ations  between  cells  are  more  sensitive  to  the  change  in  period,  suggesting 
that  the  coupling  between  cells  is  stronger. 

9.5  Accuracy  Checks  and  Computational  Time 

Figure  5-38  shows  the  results  of  a  convergence  test  applied  to  a  sinusoidal 
surface  with  d  =  0.2\,  a  =  0.  IX  and  0=0.  An  indication  of  the  accuracy 
achieved  with  an  almost  flat  surface  can  be  had  from  Fig.  5-14,  and  because  of 
the  nature  of  the  Integral  equation,  this  limiting  case  works  better  here  than  it 
did  for  E  polarization. 

Comparison  of  the  results  with  those  obtained  by  Zaki  (1969)  for  a 
sinusoidal  surface  (d  =  1.90X,  a  =  0.25X  and  0  =  0)  is  shown  in  Fig.  3-39, 
and  with  those  of  Green  (1970)  for  a  triangular  profile  (d  =  1.75X,  0^  =  20  , 

0„  =  60°,  and  0  =  12.2°)  in  Fig.  5-40.  Good  agreement  is  obtained  in  both 
instances . 
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for  H  polarization,  — o—  exact,  physical  optics. 


Fig.  5-34;  Normalized  Surface  Field  foi  a  Modified  Full-Wave  Rectified 
Surface  (d  =  0. 3X,  a  =  0. 3X.  0  »  60®  and  i  =  0. 4X)  for  H 
polarization,  — o—  exact,  —  physical  optics. 
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Fig,  5-35; 


Normalized  Surface  Field  for  a  Modified  Full-Wave  Rectified 
Surface  (d  =  0. 3X,  a  =  0. 3\,  6  =  60°  and  i  =  0,  5X)  for  H 
polarization,  — o~  exact,  —  physical  optics. 
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5-37: 
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Normalized  Surface  Field  for  a  Modified  Full-Wave 
Rectified  Surface  (d  =  0. 3X,  a  =  0. 3X,  0  =  60°  and  f  =  0. 7X) 
for  H  polarization,  — o—  exact,  —  physical  oi^ics. 
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Fig.  5-38;  Normalized  Surface  Field  for  a  Sinusoidal  Surface  (d  =  0. 2X, 
a  =  0.  IX  and  0  *=  o)  computed  using  4,  6  and  8  sampling  points, 
H  polarization 


Normalized  Surface  Field  for  a  Sinusoidal  Surface 
(d  =  1.  9X,  a  =  0. 25X  and  0  =  0),  for  H  polarization, 


(ma/  m) 
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X  A  Green 


>f— 4<  »  real  part 

ZV— O— — O  imaginary  part 


Fig.  5-40:  Surface  Field  for  a  Triangular  Surface  (d  =  1.75\,  6  =  20°, 
=  66°  and  6  =  12. 2°)  for  H  polarization 


As  in  the  case  of  E  polarization,  the  time  consumed  in  each 
computation  depends  on  several  factors  (see  section  3.5  of  Chapter  HI),  but  in 
general  the  CPU  time  is  larger  now  because  of  the  more  complicated  (modified) 
Green’s  function.  A  summary  of  the  CPU  times  for  various  combinations  of 
parameters  is  given  in  Table  V-1. 

The  distribution  of  the  CPU  time  for  a  typical  computation  where  N  =  12, 
M  =  3  and  Q  =  1  is  given  below: 

1.  Computation  of  matrix  elements  50.3  sec. 

2.  Matrix  inversion  8.4  sec. 

3.  Physical  optics  approximation  4.3  sec. 

4.  Energy  conservation  check  and  calculation  of 

diffracted  mode  amplitudes  9.0  sec. 


TABLE  V-1 

Computational  Time  for  H  Polarization 


Number  N  of 
Sampling  Points 

Number  M  of 
Terms  in  the 
Infinite  Series 

Number  Q  of 

Terms  in  the 
Periodic  Func¬ 
tion  Representation 

CPU  Time 
(sec) 

IBM/67 

4 

3 

1 

13.00 

6 

2 

10 

35.58 

6 

3 

10 

41.67 

8 

3 

1 

42.47 

8 

3 

10 

52.07 

8 

3 

12 

61.71 

8 

8 

1 

71.87 

8 

6 

1 

60.46 

10 

3 

1 

67.13 

10 

8 

1 

104.84 

12 

3 

1 

73.10 

12 

3 

10 

120.38 

16 

6 

1 

180.21 

16 

5 

1 

167.36 

20 

5 

15 

424.75 

Chapter  VT 

SCATTERED  FIELD  FOR  H  POLARIZATION 


In  this  chapter  we  consider  the  scattered  fields  scattered  energy  and 
such  associated  physical  phenomena  as  the  S  anomalies.  Numerical  data  are 
presented  to  illustrate  the  discussion,  and  comparisons  are  made  with  the 
physical  optics  approximation  as  well  as  the  results  for  E  polarization. 

6.1  Diffracted  Mode  Amplitudes 

Using  a  procedure  directly  analogous  to  that  for  E  polarization,  and 
assuming  the  incident  magnetic  field  shown  in  Eq.  (2.31),  we  can  write 


m=“Co 


A’  e 
m 


-j(8  x+X  y) 
m  m 


(6.1) 


with 

t  x'+X  y)p  -x 

^  ^  \x  -fi  Vix'nKixW  (6.2) 

I,  m  m  j 

where  X  and  ^  are  as  given  before, 
m  m 

The  dependence  on  X  and  B  is  rather  heavier  than  that  displayed  by 

m  m 

Eq.  (4.4). 

Once  the  unknown  K  (x')  has  been  found,  it  is  a  simple  matter  to  compute 

the  amplitudes  of  all  the  diffracted  modes,  both  propagating  (X  real)  and 
+  ^ 
evanescent  (Xj^  purely  imaginary),  and  the  amplitudes  of  the  propagating  waves 

for  several  different  values  of  d/X,  a/X  and  6  for  a  sinusoidal  surface  are 

given  in  Tables  VI- 1  through  VI-4.  It  is  observed  that  the  physical  optics 

approximation,  with  or  without  shadowing  included,  tends  to  underestimate 

the  amplitudes,  and  we  further  note  that  if  shadowing  is  excluded  (or  is  not 

present),  the  physical  optics  estimate  is  the  same  for  both  polarizations  (see 

Appendix  A). 
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TABLE  VI-1 


Amplitudes  of  Diffracted  Waves  for  a  Sinusoidal  Surface 
(d  =  0.2X,  a  =  0.1\). 


e 

Mag. 

1^61 

=  0° 

Phase 

0J(O) 

e 

Mag. 

hi 

=  30° 

Phase 

0  =  600 

Mag.  Phase 

hi 

Numerical 

Solution 

1 - - 

0.9830 

-12.45 

0.9910 

-23.82 

1.0640  -67.24  ! 

Physical  Optics 
w/o  shadow 

0.6425 

0 

0.7251 

0 

0.9037  0 

Ph3rsical  Optics 
w/  shadow 

0.6425 

0 

0.1550 

-65.17 

0.5440  -166.31 

TABLE  VI-2 


Amplitudes  of  Diffracted  Waves  for  a  Sinusoidal  Surface 
(d  =  l.m.  a  =  0.25X,  6  =  0). 


! 

1 

i 

j 

Mag.  Phase 

Mag. 

hi  =  Ihl 

Phase 
=  0^(0) 

Numerical 

Solution 

0.9040  -80.50 

0.3350 

52.55 

Physical 

Optics 

0.3042  0 

0.4389 

90 

TABLE  VI-3 

Amplitudes  of  Diffracted  Waves  for  a  Sinusoidal  Surface 
(d  =  0.4X,  a  =  0.2X). 


. 

0« 

0° 

0  = 

60° 

Mag. 

Phase 

Mag. 

Phase 

Kl 

0/) 

0M°> 

Numerical 

Solution 

0.9350 

-107.19 

1.0140 

-145.28 

Physical  Optics 
w/o  shadow 

0.0549 

0 

0.6450 

0 

Physical  Optics 
w/  shadow 

0.0549 

0 

-131.32 

FcroufT 


Numerical 

Solution 


Physical 

Optics 


0.9986 


0.9647 


Phase 


-0.55 


6.2  Diffracted  Enerey  and  Wood  Anomalies 


6.2.1  Discussion 

The  diffracted  energy  is  the  field  quantity  most  frequently  measured. 

It  is,  of  course,  directly  related  to  the  diffracted  mode  amplitudes,  and  the 
general  remarks  concerning  energy  and  the  Wood  anomalies  made  in  section 
4.2.1  also  apply  here.  Unlike  the  P  anomalies,  however,  the  S  anomalies 
(which  were  discovered  earlier)  have  been  more  widely  studied  and  are  more 
completely  understood.  According  to  Rayleigh's  theory,  the  occurrence  of 
the  S  anomalies  corresponds  to  the  singularity  of  the  scattered  field  at  Ray¬ 
leigh  wavelengths.  At  these  wavelengths,  a  new  spectral  order  appears  at 
grazing  angles,  and  thus  causes  a  sudden  rearrangement  in  the  diffracted 
energy  among  the  other  modes.  But  it  should  be  noted  that  there  are  S 
anomalies  other  than  iiiose  which  occur  at  the  Rayleigh  wavelengths,  and 
these  are  similar  to  the  "Brewster  angle"  type  discussed  in  section  4.2.2.  A 
comparison  with  the  P  anomalies  reveals  that  the  S  anomalies  are  generally 
much  stronger  and  can  occur  with  structures  having  relatively  shallow  grooves 
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6.2.2  Computed  Data  for  the  Diffracted  Energy 

a)  Two  radiating  modes .  The  fields  scattered  by  surfaces  of  very  small 
(electrical)  period  are  of  relatively  little  interest.  Anomalies  are  observed 
only  when  more  than  one  radiating  mode  is  present,  and,  of  course,  the  sim¬ 
plest  cases  (also  the  ones  most  Investigated)  are  tiiO?e  where  there  are  just  two 
radiating  modes. 

Figure  6-1  shows  a  plot  of  the  normalized  diffracted  energy  versus 
the  surfece  height  for  a  full-wave  rectified  profile  for  which  d  =  0.85X  and 
6  =  36°.  Since  an  equivalent  plot  for  a  sinusoidal  surface  was  obtained  by 
Zaki  (1969),  it  is  included  here  for  the  sake  of  comparison.  With  this  choice 
of  parameters,  only  the  m  =  0  (specular)  and  the  m  =  -  1  (back  scattered) 
modes  can  propagate.  The  diffracted  energy  exchange  is  continuous, 
exhibiting  maxima  and  minima  in  an  almost  standing  wave  pattern.  The  I'ate 
of  conversion  of  the  diffracted  enei^  is  more  rapid  for  the  sinusoidal  sur¬ 
face  than  with  the  full-wave  rectified  surface.  Anomalies  occur  in  both  cases, 
however,  with  the  first  appearance  being  at  a  =  0. 16X  for  the  former  surface 
and  at  a  =  0.36A.  for  the  latter.  Note  that  this  type  of  anomaly  is  different 
from  the  classical  (Rayleigh)  type  in  that  they  do  not  occur  at  the  Rayleigh 
wavelengths . 

Another  example  for  two  radiating  modes  is  shown  in  Fig.  6-2  for  a 
full-wave  rectified  profile  with  a  =  0.3\,  d  =  0. 30A,  0  =  60°  and  f  varying 
from  0  to  0.7X.  We  note  the  abrupt  exchange  of  diffracted  energy  between 
the  m  =  0  and  m=  -1  modes  near  the  Rayleigh  wavelength,  =  0.535X.. 

This  exchange  is,  however,  not  complete,  but  about  50  per  cent.  A  stronger 
anomaly  showing  complete  conversion  of  energy  is  observed  at  d  =  0.7X.. 

The  changing  division  of  diffracted  energy  between  two  radiating  modes 
as  a  function  of  the  angle  of  incidence  is  illustrated  in  Fig.  6-3  for  a  full- wave 
rectified  profile  with  d  =  0.6X  and  a  =  0.3\.  It  is  seen  that  there  is  a  dis¬ 
continuity  in  the  diffracted  energy  for  both  modes  at  0  =  41.8°  corresponding 
to  a  Rayleigh  wavelength,  and  for  incidence  angles  greater  than  41.8°  most 
of  the  energy  is  carried  by  the  m  -  -1  mode. 


s 


O— — O  sinusoidal  (Zaki) 

X - K  full-wave  rectified 


.1  .2  .3  .4  .5 

Surface  Height  (a/X) 


Fig.  6-1: 


Diffracted  Energy  for  Two  Radiating  Modes,  for 
H  polarization. 


Diffracted  Dner}-^  for  a  Modified  Full-wave  Rectilied  Surface 
(.1  O.dX,  0  =  (iO  aiul  d  =  0.  30/'\)  ,  for  II  polari/ation. 
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b)  Three  radiating  modes.  As  the  period  of  the  surface  is  increased, 
higher  order  modes  come  in,  and  the  distribution  of  diffracted  energy  among 
the  radiating  modes  becomes  more  complex.  This  is  clearly  demonstrated  in 
Fig.  6-4  for  a  sinusoidal  surface  where  d  *  1. 155\,  6  =  60°  and  a  varies 
from  0  to  0.7A.  With  this  choice  of  parameters  there  are  three  radiating  modes, 
with  the  m  *=  0  mode  radiating  in  the  specular  direction  (0^  *  60°),  the  m  =  -  J 
mode  radiating  in  the  normal  direction  (0  *  0)  and  the  m  =  -2  mode  radiating 

in  the  back  scattered  direction  2  *  “60°).  The  curves  show  no  resemblance 
to  those  for  the  case  of  two  radiating  modes.  For  a  <  0, 05\  the  m  »=  0  mode 
is  the  dominant  one,  carrying  at  least  80  per  cent  of  the  total  energy.  Once 
again,  we  see  that  the  diffracted  energy  for  each  mo<ie  is  more  sensitive  to  a 
change  in  the  surface  height  than  was  the  case  for  E-polarlzatlon,  and  a  strong 
anomaly  is  observed  for  even  the  relatively  shallow  groove  depth,  a  *  0. 3X. 

For  this  depth,  all  of  the  energy  is  carried  by  the  m  *  -2  mode  which  produces 
a  correspondingly  strong  back  scattered  field. 

The  same  kind  of  plots  were  also  obtained  using  the  physical  optics 
approximation,  either  with  or  without  shadowing,  and  they  are  shown  in  Figs. 
4-4  (the  same  for  both  polarizations)  and  6-5.  By  comparing  these  with  Fig. 
6-4,  it  is  appa'  ent  that  neither  version  of  the  physical  optics  approximation 
gives  satisfactory  agreement  with  the  exact  solution  except  perhaps  when 
a  0. 05A.. 

A  plot  showing  the  relation  between  the  diffracted  energy  and  the 
minimum  radius  of  curvature  (at  the  surface  peak  or  trough)  of  a  sinusoidal 
surface  for  which  d  =  1. 155A.  and  0  =  60°  is  shown  in  Fig.  6-6.  As  in  the 
case  of  E  polarization,  the  dominant  mode  for  ^4.5  is  the  m  =  0 

(specular)  one.  Thus,  as  far  as  the  scattered  field  is  concerned,  the  surface 

with  kp  .  >4,5  can  be  treated  as  flat. 

mm- 
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Fig.  6-4: 


Diffracted  Energy  for  a  Sinusoidal  Surface  (d  =  1. 155X  and 
0  =  60  ),  for  H  polarization. 
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yiirfact'  Amplitude  (a/ A) 


Kiii.  Dit'iraoled  I‘'noi’t>y  for  a  Sinusoidal  Surface  (d  -  1.  ir)5A  and 

'  () 

0  1)0  ).  Computed  UsiiiJi  the  Physical  Optics  Method  (witii 

shadow),  for  H  polarization. 
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6.2.3  Blazed  Gratings 

The  diffracted  energy  for  various  orders  of  radiating  mode  can  be 
computed  from  Eq.  (6.2)  and  used  as  a  basis  for  the  design  of  blazed  gratings 
(gratings  whose  profiles  are  so  shaped  as  to  maximize  the  diffracted  energy  in 
a  given  direction  or  a  given  mode). 

Blazed  gratings  find  important  applications  in  spectroscopy  where  the 
spectrum  is  observed  in  only  one  order  at  a  time,  the  energy  going  into  the  other 
order  modes  being  completely  wasted.  It  is  therefore  desirable  to  have  the 
observed  spectrum  as  luminous  as  possible,  and  one  way  of  achieving  this  is 
to  reduce  the  number  of  propagating  modes  by  using  smaller  periods.  Another 
way  is  to  minimize  the  amount  of  energy  carried  by  the  other  modes,  and 
according  to  Eq.  (6.2)  and  neglecting  the  contribution  due  to  the  perpendicular 
side  of  the  surface,  the  diffracted  energy  of  the  m^^  mode  for  a  right-angled 
triangular  profile  will  be  zero  if 


X 

m 


-)3 


m 


f'(x)  =  0 


or 


f'(x)  =  X 

m' 


m  ' 


(6.3) 


Eq.  (6.3)  is  particularly  useful  when  there  are  only  two  radiating  modes. 
However,  since  this  applies  only  to  H  polarization,  the  echelette  grating  so 
designed  is  polarization  dependent. 

6.3  Conservation  of  Energy 

One  feature  that  is  common  to  perfectly  conducting,  periodic  surfaces 
regardless  of  the  number  of  radiating  modes  involved  is  the  condition  of 
energy  conservation,  and  this  can  be  used  as  a  check  on  the  accuracy  of  the 
numerical  solution.  The  results  of  this  check  are  shown  in  Table  VI-5  for 
a  sinusoidal  surface  sustaining,  in  most  cases,  only  one  radiatkig  mode. 
Results  for  two  and  three  radiating  modes  are  shown  in  Figs.  6-1  and  6-3.  In 
all  cases  considered,  the  numerical  results  satisfy  the  condition  of  energy 
conservation,  but  the  physical  optics  approximation,  with  or  without  shadowing 
included,  does  not. 
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TABLE  VI-5 

Conservation  of  Energy  for  a  Sinusoidal  Surface 


Numerical 

Physical  Optics 

w/o  Shadow 

w/  Shadow 

d/X 

a/X 

0(°) 

Energy 

Error(o/o) 

Energy 

Error(°/o) 

Energy 

Error(°/o) 

0.2 

0.1 

0 

-1.52 

0.4728 

-58.72 

0.4728 

-58.72 

0.2 

0.1 

30 

-1.43 

0.4554 

-47.30 

0.0208 

-98.00 

0.2 

0.1 

60 

0.5075 

1.50 

0.4083 

-18.34 

0. 1479 

-70.42 

0.2 

0.01 

0 

0.9980 

-0.20 

0.9921 

-  0.79 

0.9921 

-  0.79 

0.2 

0.03 

0 

0.9987 

-0.13 

0.9308 

-  6.92 

0.9308 

-  6.92 

0.4 

0.2 

0 

0.9823 

-1.77 

0.0030 

-99.70 

0.0030 

-99.70 

0.4 

0.2 

60 

0.5042 

0.84 

0.2061 

-58.78 

0.0083 

-98.50 

1.9 

0.25 

0 

0. 9974 

-0.26 

0.4202 

-57.98 

0.4202 

-57.98 

6.4  Approximate  Techniques 

As  in  the  case  of  E  polarization,  an  attempt  has  been  made  to  use  the 
knowledge  obtained  from  the  numerical  results  to  develop  a  simple  method  for 
estimating  the  field  back  scattered  from  surfaces  at  oblique  incidence.  Were 
it  possible  to  arrive  at  a  valid  approximation  technique,  it  would  be  very  use¬ 
ful  in  such  practical  problems  as  sea  surface  or  terrain  scattering. 

Naturally  the  physical  optics  method  was  tried,  but  was  abandoned 
because  of  its  rather  poor  agreement  with  exact  results.  As  we  have  seen  in 
previous  sections,  the  scattered  field,  involving  the  conversion  of  energy  among 
radiating  modes,  is  a  very  complicated  quantity  which  offers  little  prospect 
of  simple  approximation  in  any  general  situation.  However,  in  particular  cases 
such  as  back  scattering  by  a  sinusoidal  surface  at  oblique  incidence  (Figs.  5-25 
through  5-29),  there  it^  some  hope  oi  developing  an  approximate  technique.  In 
a  procedure  directly  analogous  to  that  for  E  polarization,  the  field  back  scat¬ 
tered  by  a  sinusoidal  surface  at  oblique  incidence  was  computed  by  the  stationary 


phase  method  using  the  exact  surface  field  over  only  a  restricted  portion  of 
the  illuminated  region  of  the  surface.  The  results  are  shown  in  Table  VI-6. 


TABLE  VI-6 

Comparison  of  the  Energy  Backscattered  by  a 
Sinusoidal  Surface:  d  *  1. 155>,  6  -  60^,  Computed 
Using  the  Exact  and  Stationary  Phase  Methods. 


Surface 

Amplitude  (a/\) 

Normalized  Back  Scattered  Energy 

Exact  Method 

Stationary 

Phase  Method 

0.3 

0.98 

1.06 

0.4 

0.64 

0.69 

0.5 

0.75 

0.70 

0.6 

0.60 

0.64 

0.7 

0 

0.16 

Although  the  agreement  with  the  exact  values  is  excellent,  the 
approximate  method  still  requires  a  knowledge  of  the  true  surface  field  over 
a  small  region  of  the  surface,  and  in  order  t'  avoid  solving  the  integral  equa¬ 
tion  numerically,  it  is  necessary  to  find  some  other  means  for  estimating  this 
field.  As  regards  the  phase,  the  physical  optics  value  turns  out  to  be  adequate; 
unfortunately,  the  mod'll  us  is  a  much  more  complicated  quantity  and,  unlike 
E  polarization,  is  not  dependent  on  the  local  surface  properties.  Many  different 
methods  using  combinations  of  ray  optics  and  surface  wave  theory  have  been 
tried.  For  example,  one  method  considered  the  Interaction  of  two  rays:  the 
incident  ray  and  the  ray  reflected  from  the  opposite  face  of  the  surface;  another 
method  considered  forward  and  backward  traveling  v/aves  whose  amplitudes 
we  sought  to  choose  to  fit  the  exact  data.  However,  neither  of  these  methods 
alone  (or  even  in  combination)  was  satisfactory,  though  it  was  noted  that  a 
rough  estimate  for  the  modulus  was  provided  by  the  imaginary  part  of  the 
physical  optics  approximation.  Though  there  is  no  physical  basis  for  this 
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whatsoever,  the  corresponding  estimates  of  the  back  scattered  field  were 
computed,  and  a  comparison  with  the  exact  values  is  given  in  Table  VI-7  . 


TABLE  VI-7 

Comparison  of  the  Current  Modulus  over  the  Illuminated  Region 
Of  a  Sinusoidal  Surface,  Computed  Using  the  Exact  Method  and 
The  Imaginary  Part  of  the  Physical  Optics  Approximation. 

(d  =  1.155X,  e  =  60°). 


Normalized  Current  Modulus 

Location 

Approximation  Using 

Surface 

of  Point, 

Exact 

The  Imaginary  Part  of 

Amplitude  (a/X) 

X  Coordinate  (X) 

Method 

The  Phys.  Optics  Values 

0.915 

1.75 

2.00 

0.3 

1.010 

2.05 

1.92 

1.109 

2.30 

1.84 

0.4 

0.938 

1.84 

1.96 

:  1.082 

1.99 

2.00 

’  0.915 

2.00 

1.97 

0.5 

1.010 

1.42 

1.90 

i  1.109 

1.84 

1.96 

0.915 

1.85 

1.9C  ! 

0.6 

,  1.010 

1.17 

1.70  1 

1.109 

1.84 

1.70 

0.7 

j  1.010 

1.04 

1.43 

i  1.109 

0.90 

1.37 

Chapter  VII 

RAYLEIGH'S  METHOD 

As  mentioned  in  Chapter  IV  in  connection  with  Wood  anomalies, 
Rayleigh's  work  was  one  of  the  early  attempts  to  theoretically  investigate 
scattering  (electromagnetic  or  acoustic)  by  a  periodic  surface.  His  approach, 
now  termed  "Rayleigh's  method",  is  indeed  an  important  one,  and  has  been 
used,  with  or  without  refinement,  by  many  subsequent  investigators.  We 
shall  now  examine  his  method  in  some  detail  and  compare  exact  data  with 
data  obtained  using  his  assumption. 

7. 1  History 

In  11)1)7  Lord  Rayleigh  studied  the  scattering  of  sound  waves  by  corru¬ 
gated  surfaces  and  sought  to  solve  the  wave  equation  in  combination  with  the 
appropriate  boundary  condition.  His  method  is  an  intuitive  one  that  has  been 
used  by  many  mvestigators  and  is  based  on  the  following  two  assumptions: 
a)  the  scattered  field  may  be  represented  by  a  linear  combination  of  discrete 
plane  waves,  either  propagating  or  attenuated  from  the  surface,  and  b)  this 
representation  holds  everywhere  above  and  on  the  boundary.  A  simp!“  ^nd 
straightforward  description  of  the  Rayleigh  method  for  a  periodic  boundary  is 
given  by  Beckmann  and  Spizzichino  (1963).  The  second  assumption,  which  is 
often  termed  the  "Rayleigh  assumption,  "  has  been  a  controversial  subject  for 
many  years . 

As  first  observed  by  Llppmann  (1953),  a  choice  of  scattered  field 
consisting  only  of  outgoing  waves  is  not  the  most  general,  and  may  not  com¬ 
pletely  represent  the  field  within  the  grooves;  at  any  point  in  this  region,  one 
would  expect  the  field  to  Include  waves  propagating  (or  exponentially  damped) 
In  both  directions,  since  there  are  currents  on  those  parts  of  the  grating  above 
and  below  the  given  point.  Meecham  (1956b)  studied  reflection  from  irregular 
surfaces  with  a  method  based  on  Fourier  transforms,  and  obtained  numerical 
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results  for  the  distribution  of  energy  reflected  from  a  sinusoidal  surface.  By 
comparing  these  with  experimental  results  and  also  with  the  ones  obtained  using 
Rayleigh's  method,  he  concluded  that  errors  were  of  the  same  order  as  those 
in  the  physical  optics  approximation. 

Heaps  (1957)  attacked  the  validity  of  the  "Rayleigh  assumption"  in 
another  way.  He  studied  the  reflection  of  plane  acoustic  waves  at  a  sinusoidal 
surface  of  zero  pressure  and  investigated  the  least  possible  value  of  the  sur¬ 
face  pressure  consistent  with  the  assumption  that  all  the  reflected  radiation 
is  in  the  form  of  undamped  plane  waves.  After  comparing  his  results  with 
the  experimental  data  of  LaCasce  and  Tamarkln  (1956),  he  concluded  that 
"if  all  reflected  energy  has  the  form  of  undamped  plane  waves  then  the  surface  is 
necessarily  sound  absorbing  and  of  pressure  significantly  different  from  zero. 
Thus,  in  the  neighborhood  of  the  corrugated  surface  of  zero  pressure,  it  is 
necessary  to  take  into  account  other  forms  of  radiation  and  such  forms  play  a 
significant  part  in  satisfying  the  boundary  condition.  "  Murphy  and  Lord  (1964) 
have  also  criticized  tlie  Rayleigh  assumption.  By  making  a  direct  comparison 
of  the  results  obtained  from  an  exact  representation  of  the  field  with  those  based 
on  Rayleigh's  formulation,,  they  concluded  that  the  latter  is  inadequate. 

On  the  defense  side,  i  ^re  aie  people  who  support  Rayleigh's 
assumption.  Marsh  (1963)  has  generalized  Rayleigh's  method,  and  considered 
the  scattering  of  acoustic  waves  from  a  sinusoidal  surface  (y  *  a  cos  Kx)  with 
a  Dirlchlet  boundary  condition.  He  showed  that  the  boundary  condition  can  be 
satisfied  using  only  outgoing  waves.  Numer’ical  data  were  obtained  for  the 
amplitudes  of  diffraction  which,  upon  substitution  into  the  energy  relation,  in¬ 
dicated  conservation  of  energy.  However,  it  must  be  noted  that  the  values  of 
Ka  used  were  relatively  small,  and  as  shown  later,  the  Rayleigh  assump¬ 
tion  is  indeed  valid  In  this  case.  Additional  support  of  Rayleigh's  assumption 
comes  from  Millar  (1963),  who  argued  from  the  viewpoint  of  analytic  continu¬ 
ation.  He  i'emarked  that  if  the  boundary  values  are  analytic,  "the  exterior 
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field  may  be  continued  into  the  grating  from  above,  and  that  the  exterior  plane 
wave  representation  will  converge  everywhere  above  a  plane  parallel  to  y  =  0. 
Likewise  uie  interior  field  may  be  continued  into  the  upper  half-space  from 
belc  and  the  interior  plane  wave  representation  will  be  convergent  below  a 
plane  parallel  to  y  =*  0,  If  the  two  half-spaces  of  convergence  so  defined  con¬ 
tain  a  common  layer  which,  in  turn,  completely  contains  the  reflecting  surface 
y  ==  a  cos  Kx,  then  the  two  representations,  each  outgoing  or  attenuated  in  its  own 
half-space,  may  be  matched  on  the  reflecting  surface  where  both  ax'e  valid.  " 
This  remark  seems  very  reasonable. 

Possibly  the  most  significant  step  toward  the  resolution  of  the  contro¬ 
versy  (with  Dirichlet  condition)  was  taken  by  Petit  and  Cadilhac  (1966).  They 
demonstrated  (using  analytic  continuation)  that  for  the  profile  y  *  a  cos  Kx 
(-00  <  X  <  oo)  Rayleigh's  assumption  could  not  be  valid  if  Ka  >  0. 448;  for 
Ka  <  0.448,  the  question  remained  open.  Recently  Millar  (1969)  used  a  tech¬ 
nique  developed  to  locate  singularities  of  solutions  to  the  Helmholtz  equation. 

By  applying  this  to  the  sinusoidal  profile,  he  was  able  to  show  that  Rayleigh's 
assumption  is  indeed  valid  if  Ka  <  0. 448  and  is  not  valid  if  Ka  >  0. 448. 

In  order  to  shed  some  light  on  this  remark,  let  us  compute  the  scattered 
field  using  the  exact  representation  and  under  the  Rayleigh  assumption. 

7 .  2  Numerical  Results  by  the  Exact  and  Rayleitdi  Methods 
For  E  polarization  we  have,  from  Eq.  (2.22), 

E^(x,y)==^  \  P  \/l  f '(x)^  K  (x)dx'  .  (7.1) 

7.  4^1  z 

Jo 

Using  the  exact  representation  we  can  expand  the  field  in  terms  of  incoming  and 
outgoing  waves,  thus 
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E  (x,y)  =  ✓  .  B  e  C  e 
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(7.2) 


From  Eqs.  (7. 1)  and  (7.2)  we  get 
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\/l  +  f(x)^  K  (x)dx  ,  (7,  3) 


-X^y)  I - 2 

Vl+f*(x)  K  (x)dx 
z 


(7.4) 


where  S  and  S  represent  the  portions  of  the  surface  above  and  below  a  plane 
through  the  point  P(x,  y)  (see  Fig.  7-1).  Of  course,  when  y  is  above  the  grating, 
i.e.  y  >  max  f  (x)  ,  then 

B  *  A  and  C  =  0  , 

mm  m 

where  A  is  defined  by  Eq.  (4. 4) 
m 

On  the  other  hand  using  Rayleigh’s  assumption,  we  can  expand  the  field 
ever3ovhere  in  terms  of  outgoing  waves  alone,  viz. 


00 
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If  we  now  substitute 
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into  Eq.  (7. 1),  we  get 
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Ig  m=-OD  m 


K  (rjdx. 
z 


(7.6) 


In  order  to  make  Eq.  (7. 5)  equal  to  Eq.  (7. 6)  everywhere  on  the  surface,  w’e 
must  remove  the  modulus  sign  in  the  kernel.  Thus,  equating  Eqs.  {?.  5)  and 
(7. 6),  we  obtain 
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Vl+fCx)"  K^(x)dx  .  (7.8) 


Intuitively  we  expect  the  discrepancy  between  Eqs.  (7. 2)  and  (7. 5)  to  be  greatest 
when  S  =0,  that  is  when  the  point  P(x,  y)  is  In  the  surface  trough.  Thus,  let¬ 
ting  S  =  0,  we  get  from  Eq.  (7. 2) 
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and  from  Eq.  (7.5) 
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where 


B  » 
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Using  Eqi .  {7. 9)  and  (7. 10>,  numerical  results  for  the  scattered  field  were 
computed  for  a  sinusoidal  and  an  inverted  full-wave  rectified  surface,  and 
these  are  shown  in  Table  Vn-l. 


TABLE  VH-l 


E-polarization.  Comparison  of  Scattered  Fields  at  the  Surface 
Trough  Obtained  by  the  Exact  and  Rayleigh  Mefiiods. 


Type  of 
Surface 


Sinusoidal 

Inverted 

Full-wave 

Rectified 


Scattered  Fj 

[eld 

Ka 

Exact  Representation 

Raylei^'s  Method 

0. 314 

0. 942 
1.884 

2.27440  -jO.  22657 
1.96310  -jO.  48310 

1. 86610  -  jO.  52445 

2. 27420  -  jO.  22664 
1.95450  ~j0. 48166 

1. 75170  -  jO.  54162 

0.314 

0.942 

2.42170  -jO.  18443 

2. 22550  -  jO.  49810 

2.42180  -jO.  17913 
2.22720  -jO.  46486 

Table  VII-l  shows  that,  for  sufficiently  small  values  of  Ka,  the 
scattered  fields  based  on  the  exact  and  Rayleigh  methods  agree  extreme./ 
well  for  a  sinusoidal  surface.  This  is  consistent  with  Millar's  remark.  In 
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the  case  oi  the  inverted  full-wave  rectified  surface,  agreement  is  rtot  so  good 
for  the  same  values  oi  Ka,  This  is  probably  due  to  die  fact  that  fee  ir-re'-ted 
full-wave  rectified  profile  is  not  an  analytic  function  as  is  requii  ad  hr  iiiillar's 
analysis. 

Turning  now  to  H-polarization,  we  have 
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Using  the  exact  representation  we  get 
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and  when 


y_>  max[f{x)]  , 


then 


B'  =A» 
m  m 


and 


C  =0 
m 


where  is  defined  by  Eq,  (6-2). 
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On  tbs  other  hand,  using  l?a5ietgti*8  assumption,  we  !mve 


mm 

H  (x,y)=  y  ,  B*  e 
z  ■'  A—*  m 
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(7. 15) 


If  we  now  substitute  fee  expression 
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into  Eq.  {7, 11)  we  get 
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Again,  in  order  to  make  Eq,  (7. 15)  equal  to  Eq.  (7. 16)  everywhere  on  the 
surface,  we  require  the  modulus  sign  in  the  kernel  to  be  removed.  Thus, 
equating  Eqs.  {7, 15)  and  (7. 16), 
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and  if  v/e  examine  the  field  at  the  surface  trou^  (so  that  S  ^  0)^  then 
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(7. 20) 


Numerical  values  computed  using  Eqs.  (7. 19)  and  (7. 20)  for  a  sinusoidal 
surface  and  an  inverted  full-wave  rectified  surface  with  d  =  0. 2\  0*0  and 
different  values  of  Ka  are  shown  in  Table  Vn-2. 


TABLE  Vn-2 

H-polarization.  Comparison  of  Scattered  Fields  at  Surface 


Trough  Obtained 

by  the  Ex,>ct  and  Rayleigh  Methods 

Type  of 

Scattered  Field 

Surface 

Ka 

Exact  Representation 

Rayleigh's  Method 

Sinusoidal 

0.314 

0.157 

0. 74072  -  jO.  63363 
0.87126  -jO.  31756 

0.  75254  -  jO.  62012 
0.87870  - jO.  31196 

Inverted 

Full-wave 

Rectified 

0.314 

0.942 

0. 79608  -  jO,  34601 
-0.07901  -jO.  90223 

0. 92741  -  JO. 34870 
1.07230  -  Jl.  00060 
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It  is  interesting  to  note  that  the  critericHi  under  ^*diich  fee  Raylei^ 
assumpticm  holds  seems  different  for  H-polari2ati<m,  wife  all  data  su^esting 
that  Ka  must  be  smaller  than  the  value  0. 448  appropriate  to  E-poIarization. 
However,  the  exact  figure  for  H-polarization  is  unknown.  Again  the  agree- 
mept  between  fee  exact  representation  method  and  the  Rayleigh  me&od  be¬ 
comes  worse  if  the  periodic  surface  Is  not  an  analytic  function  (sueb  as  an 
inverted  full-wave  rectified  profile). 

From  the  above  computed  data  for  E  and  H  polarisations,  we  there¬ 
fore  conclude  that  Raylei^’s  mefeod  can  be  used,  with  sufficient  accuracy, 
for  estin^cting  fee  field  scattered  by  a  periodic  surface  provided  feat  it  is 
analytic  and  has  a  email  value  of  Ka  (say  less  then  0. 5).  If  the  surface  is  not 
an  analytic  curve  {in  practice,  a  grating  profile  cannot  be  shaped  to  coincide 
with  an  analytic  curve),  Rayleigh’s  mefeod  providesi  us  better  results  if  fee 
incident  field  is  E-polarized. 


Chapter  VIE 

APPLICATION  TO  ROUGH  SURFACE  SCATTERING 
8. 1  Discussion 

Although  the  main  purpose  of  this  work  is  concerned  with  scattering 
by  a  periodic  surface  of  arbitrary  shape,  it  is  appropriate  to  consider  the 
application  of  our  results  to  scattering  from  rough  surfaces. 

In  the  past  decade  or  so,  a  great  number  of  papers  have  been  published 
dealing  with  this  type  of  scattering.  A  few  of  the  more  important  refercaoss 
are  Rice  (1951),  Hoffman  (1955),  Beckmann  (1963)  and  Fung  and  Ch  ui  (1939), 
and  if  there  is  one  feature  that  almost  ail  such  theoretical  treatments  have 
in  common,  it  is  their  reliance  on  the  physical  optics  approximation  which,  in 
effect,  reduces  the  determination  of  the  scattered  field  to  quadratures.  Of 
course,  there  exist  many  differences  of  detail,  both  conceptual  and  analytical, 
depending,  for  example,  on  whether  the  perturbed  surface  is  defined  specifically 
or  statistically,  the  nature  of  the  statistical  variations  assumed,  the  stage  at 
which  averaging  is  performed,  the  method  of  evaluating  the  remaining  surface 
integral,  whether  shadowing  is  taken  into  account,  etc.  However,  the  results 
obtained  are  still  limited  by  the  physical  optics  approximation,  which  takes  no 
account  of  any  multiple  scattering,  is  invalid  for  surfaces  of  curvature  small 
compared  to  the  wavelength  (small  scale  roughness),  and  provides  a  false  esti¬ 
mate  of  shadowing. 

If  the  surface  is  illuminated  near  normal  incidence,  these  deficiencies 
are  of  little  consequence,  primarily  because  the  base  surface  provides  a  strong 
return.  As  we  depart  from  normal  incidence,  the  surface  concentrates  the 
return  in  the  forward  (specular)  direction,  where  the  estimated  scattering 
from  the  perturbed  surface  may  remain  valid,  but  in  other  directions  (e.g. 
back  scattering)  all  of  the  observed  return  is  a  consequence  of  roughness,  and 
any  perturbation  based  on  physical  optics  is  of  doubtful  accuracy. 
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In  order  to  treat  adequately  the  case  of  oblique  incidence  on  a  rough 
surface,  it  is  desirable  to  postulate  a  base  surface  which  will  itself  generate 
a  back  scattered  return  at  these  angles;  to  devise  a  means  for  accurately 
estimating  this  scattering;  and  then  to  handle  by  a  peiturbation  technique  such 
small  scale  roughness  as  it  may  be  desirable  to  includa.  In  order  to  simulate 
actual  rough  surfaces,  the  small  scale  roughness  will  l»e  assumed  statistical 
in  character.  From  the  viev;poiat  of  optical  gratings,  the  effect  of  this  kind  of 
roughness  is  equivalent  to  degrading  a  diffraction  grating  '^dth  small  random 
errors.  For  this  reastm,  the  results  obtained  here  are  also  applicable  to  the 
study  of  grating  imperfections  which  may  cause  phenomena  such  as  "ghosts" 
and  wave-front  aberrations. 


8. 2,  Perturbation  of  the  Surface  Field 

Consider  a  plane  wave  of  either  polarization  incident  in  a  direction 
making  an  angle  d  with  the  positive  y  axis,  as  shown  in  Fig.  8-1.  We  repre¬ 
sent  the  rough  surface  by  a  small  scale  ioughness  superimposed  on  a  sinusoidal 
base  simulating  large  scale  roughness.  The  surface  height  is  then 

y  =  fix)  =  acos^^“x^  +  5(x)  ,  -a)<x<oo  (8.1) 

with  a»|5(x)|. 

Furthermore  we  assume  that  6(x)  is  a  random  variable  obeying  a 
Gausslaii  distribution  with  zero  mean,  such  that 


g 


exp  (- 


(8  2) 


<fi(x)>=  0 

where  o  is  the  standard  deviation,  and  the  bracket  sign  <  >  represents  a  statis¬ 
tical  average  process. 
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In  accordance  with  the  boundary  conditions,  the  wave  induces  on  the 
surface  a  current  density  F(r)  which  is  the  actual  source  oi  the  scattered 
field.  Application  of  Green’s  second  identity  and  the  appropriate  boundarj' 
conditions  then  yields  the  following  expressions  for  the  scattered  field: 
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for  £ -polarization  and 


aOO 

H«(,, ,  .JaSli  \  ,JWxcos0+ JstaW 

Z  CJ/LI  1  p 

J-oo 


kq(r) 

u/u 


for  H-polarlzatlon 


where 


r® 

l  cos  0- sill  0|  F^(x)dx 

J-oo 


(8.6) 


Since  the  perturbation  is,  by  assumption,  small,  the  perturbation  on 
the  periodic  base  is  equivalent  to  a  perturbation  on  the  surface  current.  More¬ 
over,  by  using  the  fact  that  the  moduluo  of  the  surface  field  is  less  sensitive 
to  a  small  change  in  the  surface  height  than  the  {^ase,  we  can  assume  that  the 
perturbation  on  the  surface  field  produces  only  a  small  change  in  the  phase  of 
the  unperturbed  surface  field.  This  allows  us  to  replace  the  rough  surface  by 
its  periodic  base.  Let 


y=sf(x)»aco8(^-^x^ 


and 


for  E-polarlzation, 


F  (x)»K  (x)e 
z  z 


jWx) 


(8. 7) 


F^(x) «  K^(x)e' 


jka(x) 


(8.8) 


for  H-polarization,  where  K  (x)  and  K  (x)  are  the, surface  fields  for  the  un.ier- 

z  t 


turbed  surface. 
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Substituting  Eq,  (8.  ?)  into  Eq.  (8. 5)  we  get 


e\,  0)  -  q(r)  \ 

J-oo 


dx  (8.9) 


The  intensity  is  then 


E®(r,0)p  =  E®(r,0)*E®%,0) 

Z  I  z  z 


=iq(r)f 


A®  A® 

4'-00  J-oc 


jk[(x,-X2)cos0+(yj-y2)sin0]  x 


(l  +  f^  jl+f(xj^  Kjx,  )Kjx„)  e 


jk[^Xj)-6(x2^ 


2  ‘*2  1  z'  2 


dXjdXg, 


(8. 10) 


which  has  the  average  value 


p®  ^®  jic[(x^-X2)eo8  0+(yj-y2)sinj^^ 


E%.0)i^^^q(r)f\  \ 


•00  ^-00 


K  (x  )K(x  )<e 
Z  X  z  ^ 


jk  |^Xj)-6(x2^ 


>dxidx2 


(8.11) 


From  our  knowledge  of  ataiistics  it  is  known  that  if  fi(x. )  and  fi(x  )  are  two 
random  variables  of  normal  distribution  with  zero  mean,  variance  a  and 
correlation  coefficient  cii),  their  joint  probability  function  is  (see,  for  example, 
Middleton,  1980) 


f|6(Xj^),(5(x2)| 


2  2 
2jra  a -c) 


f”  6(Xj^)^-2cfi(x^)fi(x2)  +  6(x2)^ 


2o^{l-c^) 


(8.12) 


t 


The  factor  <  >  in  Eq.  (8.  U)  can  now  be  evaliwted  amd  Is  (Beckmann  and 
Spizzichincj  iS63) 

jkfex,)-a[x  >j  r  9  9 

<e  ^  '^■'>*exp|-kV{l-c)j'  . 

and  when  this  is  subsdtutsd  into  Eq.  (8. 11),  we  obtain 


(8. 13) 


<  E^{r 
z 


i^oo  ^0^ 

.0)f>=lq!rtf\  I  > 

v-ooj-co 


jk^Xj-X2)cos  0+(y^-y,)sin^ 


•JK  (x,  )K  (x_)exp  1  -  cT|  dx  dx  (8. 14) 

ZiZZ  I-  1  d 


For  convenience,  the  normalized  correlation  coefficient  will  be  used 
c(t)  =*  exp(-T^/T^) 


(8. 15) 


where  r=x^-x^  (assuming  stationary  process)  and  T  is  the  correlation  distance. 


Using  the  expansion 


,  r52,,  n  -nrV 

exp  [-k  cT  (1-c^  -  e  e 


a  little  mathematical  manipulation  is  sufficient  to  convert  Eq.  (8. 14)  Into 


If®  2  if® 

<  I  E®(  r.  0)|  v(r)(T/ka)^  \  u(x,  0)dx  +  2v(r)  \  xhx.  ^)it  e 

J-co  J-ao 


-x^/T^  -2 


+  2v(r) 


S.oo 

u(x, 

*  -on 


2  -x^/T^  2 

0)x  s  '  dx  +. 


(8. 16) 


5 


v/here 


v(r)  =Jq(r)r  (lor/T)  , 

lyx) . 

To  evaluate  Eq.  (8. 16)  numerically,  we  must  reduce  the  limit  of 

integration  from  infinity  to  a  finite  value.  This  can  be  done  by  noting  that 

K  (x)  is  a  periodic  function.  Thus,  using  Eq.  (2. 21), 
z 


X",  0)P  >  *  v(  r)(T/kcr)^ 


+  2v(r) 


I  Jo  “=-<» 


jkmd(cos  0  -  sin  J 


\  u(x,0)  ^  /T  jdxj 


+  2v<r) 


5  ?  2  9 

U(x,0)  Idx' 

0  ^ 


(8.17) 


which  can  be  evaluated  by  truncating  the  infinite  series.  On  examining  Eq. 

(8. 17),  it  is  observed  that  for  small  values  of  a  the  first  term  dominates  and 
represents  contributions  coming  mainly  from  the  periodic  base,  while  the 
second  and  higher  order  terms  are  the  contributions  from  the  perturbations- 
Numerical  results  showing  the  scattered  field  pattern  for  a  sinusoidal 
profile  (d  =  1. 155A,  a  =»  0. 3X  and  0  =  60^)  are  shown  in  Figs.  8. 2  through  8. 4 
with  a  =  0. 05\  and  T  varying  from  0. 5d  to  3d.  It  Is  seen  that  the  diffracted 
energy  spectrum  has  broadened  slightly  into  a  shape  like  a  pencil  beam  around 


Fig.  8-6: 


Scattered  Intensity  Pattern  (T  =  d),  for  H  polarization. 
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the  discrete  (unperturbed)  spectrum,  bi  all  cases  considered,  there  is  a 
decrease  of  diffracted  energy  in  the  specular  and  normal  directions  when  the 
surface  is  perturbed.  There  is  also  some  evidence  of  surface  wave  radiation 
for  small  0.  In  the  back  scattered  direction,  however,  the  energy  is  relatively 
unaffected. 

In  an  analogous  way,  using  Eq.  (8. 6)  we  get  for  H-polarization 


H®(r,0)p>»v'(r)(T/lOT)" 


pd  00 

I  \  u*(x,  0) 

1  m=-<x>  ^ 


jkmd(cos  0  -  sin  0)' 


r  ^  2  2 

■2v'(r)  \  u';*.0)  j-dx 

Jo 

^  2  2  i2 

■2v'(r)  \  nkmd(oo30-sto0),^^_^^,2^-(x+md)  /T 

Jo 


+  .  .  .  . 


(8.18) 


where 

v'(r)  •"  (k/u)^v(r) 

u'(x,  0)  =  e^^^^°°®^^^^^”^^-|^f'(x)cos0-sin0j'  K^(x)  , 

in  which  the  first  term  represents  the  dominant  contribution  from  the  periodic 
base  surface,  and  the  second  and  higher  order  terms  show  the  effect  of  the 
small  scale  roughness.  Numerical  results  for  the  scattered  field  pattern  for 
a  sinunoidal  surface  (d  =  1, 155X,  a  =  0.4\  and  0  =  60°)  are  presented  in  Figs. 
8-5  and  8-6,  where  0  =  0. 05\  and  T  varies  from  0. 5d  to  d.  It  is  seen  that 
they  bear  the  same  features  as  those  for  E-polarizatlon;  however,  we  observe 
now  there  is  a  decrease  of  diffracted  energy  in  the  back  scattered  direction. 
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In  the  above  study  we  have  assumed  an  unlimited  illumination  area  viiich 
means  that  the  entire  infinite  rough  surface  is  exposed  to  the  incident  wave 
(whether  there  is  shadowing  involved  is  another  matter).  However,  in  a  mors 
practical  situation,  it  is  likely  that  not  all  of  the  rough  surface  will  be  illuminated 
(limiting  illumination  area)  because  of  the  finite  aperture  of  the  illuminating  source. 
Thus  let  us  assume  the  incident  wave  to  be  a  plane  wave  vdiose  amplitude  is  given 
by  a  Gaussian  distribution  as 


E‘(r)  =  I(x)e'W’‘ 

Z  — 

for  E -polarization,  and 

z 

for  K- polarization,  where 


(8. 19) 


(8.20) 


I(x)  =  exp(-xV2p^)  (8.21) 

and  p  is  the  standard  deviation. 

Rigorously  speaking,  in  order  to  find  the  current  distributions  on  the 
rough  surface  due  to  the  finite  source  illuminations  Eq.  (8. 19)  or  (8. 20), 
it  is  necessary  to  solve  two  other  problems  whose  integral  equations  are 


for  H-polarlzatlon. 
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But  if  the  amplitude  of  the  Incident  wave  doea  not  decrease  very  rapidly  with 
respect  to  x,  the  current  distribution  for  the  finite  source  illumination  can  be 
approximated  by  the  current  distribution  for  the  non-limiting  source  illumination 
after  multiplication  by  the  factor  given  in  Eq.  (8. 21).  ITi'-is  from  Eq.  (8. 16) 


ico 


E®(r,  0)|^ > *  v{r)(T/ko)^ 


u(x,  0)d:{ 


'-00 

"kOO 


+  2v<r) 


^(x,  0)xe  ^  dx 


/-oo 


+  2v(r) 


i 


(8.24) 


where  u(x,  0)  =«  u(x,  0)I(x). 

Eq.  (8, 24)  can  be  computed  numerically  by  considering  only  that  part 
of  the  surface  where  the  amplitude  of  the  Incident  wave  is  significant.  The 
numerical  results  for  a  sinusoidal  surface  (d  ®  1. 155X,  a  *  0.  3X,  and  0  *  60°) 
with  p  =  3d  and  the  limit  of  integration  varying  from  -lOd  to  +10d  are  shown 
in  Fig.  8r7.  It  is  seen  that  the  results  are  very  similar  to  those  obtained 
previously  for  an  infinite  illumination  area. 

in  an  analogous  way,  we  obtain  for  H-polarization 


r 

J 
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2  X  V  , 

®  dx 


(8.25) 


where  'u’(x,  0)  *  u'(x,  0)I(x). 


The  numerical  results  are  shown  in  Fig.  8. 8  (using  the  same  parameters 
as  for  E“polarizavion).  Again  it  is  seen  that  the  scattered  field  intensity  pattern 
has  the  same  characteristics  as  the  pattern  for  an  infinite  illumination  source. 
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Chapter  IX 
CONCLUSION 

We  have  here  considered  the  problem  of  a  plane  electromagnetic  wave 
incident  on  a  perfectly  conducting,  two  dimensional  i^ariodic  surface,  and  have 
developed  numerical  procedures  for  the  direct  digital  solution  of  the  Integral 
equations  for  the  surface  'Ab.  By  using  special  summation  techniques  fol¬ 
lowed  by  the  subtraction  of  the  dc  term  to  Improve  the  convergence  of  the 
series  for  the  modified  Green's  function,  a  relatively  efficient  procedure  has 
been  arrived  at,  and  this  has  been  programmed  for  a  computer. 

For  both  polarizations  (E  and  H),  data  have  been  obtained  for  the  surface 
fields  on  a  variety  of  profiles  representing  sinusoidal,  full-wave  rectified, 
inverted  full-wave  rectified  and  triangular  surfaces  having  different  periods, 
amplitudes  and  angles  of  incidence.  In  each  case  the  physical  optics  approxi¬ 
mation  is  included  for  comparison  and  to  aid  in  developing  a  simpler  method 
for  predicting  scattering  from  a  periodic  surface.  It  was  found  that  the  polari¬ 
zation  has  a  marked  effect  on  the  surface  field,  and  the  current  distribution  is 
strongly  dependent  on  the  number  of  radiating  modes.  In  general  the  physical 
optics  method  fails  to  give  an  accurate  prediction. 

Knowing  the  surface  field,  the  amplitudes  of  the  diffracted  waves  in  the 
discrete  angular  spectrum  representation  of  the  scattered  field  can  be  computed, 
and  this  has  been  done  using  the  exact  surface  fields  as  well  as  the  physical 
optics  estimates  with  shadowing  either  included  or  Ignored.  Here  again  the 
physical  optics  predictions  are  deficient,  and  whereas  the  results  derived  from 
the  numerical  program  satisfy  the  conservation  of  energj'  law,  the  physical 
optics  values  do  not. 

The  numerical  results  have  been  used  to  explore  the  physical  phenomena 
associated  with  the  diffracted  energy  of  the  radiating  modes  and  known  as  Wood 
anomalies.  Both  the  S  and  P  anomalies  were  observed.  It  was  found  that  there 


are  anomalies  which  do  not  occur  at  the  Rayleigh  wavelengths,  and  the  S-type 
anomalies  are  generally  stronger  than  the  P-type  and  can  occur  on  structures 
having  relatively  shallow  grooves. 

Consideration  has  also  been  given  to  the  Rayleigh  assumption.  Numerical 
data  for  the  scattered  field  were  obtained  for  surfaces  of  different  profiles  using 
the  exact  and  Rayleigh  methods.  The  results  are  consistent  with  Millar’s  re¬ 
mark,  and  indicate  that  the  criteria  under  which  the  Rayleigh  assumption  holds 
are  different  for  different  polarizations. 

Lastly  the  knowledge  gained  from  the  study  of  periodic  surfaces  was 
applied  to  a  study  of  rough-surface  scattering.  The  rough  surface  was  con¬ 
sidered  as  a  small  scale  roughness  superimposed  on  a  sinusoidal  periodic  base 
(large  scale  roughness).,  the  small  scale  roughness  having  a  Gaussian  distribu¬ 
tion.  The  numerical  results  obtained  show  that  the  main  effect  of  a  small  per¬ 
turbation  is  to  broaden  slightly  the  discrete  (unperturbed)  spectrum. 
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Appendix  A 


PHYSICAL  OPTICS  APPROXIMATION 

Although  the  main  purpose  of  this  work  is  the  development  of  a  numencal 
technique  for  the  determination  of  exact  values  for  the  surface  and,  hence,  scat¬ 
tered  fields,  we  have  found  it  desirable  to  compare  the  results  obtained  with 
those  provided  by  the  physical  optics  aHproximation.  This  approximation  is 
a  ph5r3ically-based  one  which  postulates  an  explicit  form  for  the  suriace  field 
arrived  at  by  assuming  each  element  of  the  surface  to  bear  that  current  which 
it  would  carry  were  it  part  of  the  local  tangent  plane.  The  calculation  of  the 
scattered  field  is  then  reduced  to  quadratures. 

In  many  instances,  however,  and  a  periodic  sheet  is  one,  an  anals^tical 
evaluation  of  the  physical  optics  integral  is  a  difficult  procedure,  particularly 
in  such  cases  where  part  of  the  surface  is  shadowed,  and  some  of  the  short¬ 
comings  of  pitysical  optics  estimates  in  general  certainly  arise  from  sloppiness 
in  the  evaluation  of  the  integral.  In  still  other  cases,  the  physical  optics  result 
proves  to  be  more  accurate  if  shadowing  is  ignored  (see,  for  example,  Adachi, 
1965)  and  it  is  of  interest  to  observe  that  for  a  sinusoidal  sheet  it  is  then  pos¬ 
sible  to  produce  an  exact  evaluation  of  the  integi-al.  Needless  to  say,  however, 
the  resrJting  scattered  field  is  still  subject  to  the  unknown  errors  inherent  in 
the  use  of  the  physical  optics  approximation,  and  to  the  neglect  (if  present)  of 
all  shadowing  effects. 

The  procedure  is  directly  analogous  to  that  given  by  Senior  (1959)  for  the 
particular  case  of  an  H  oolarized  plane  wave  at  normal  incidence  on  a  sinusoludl 
sheet  and  consists  of  three  steps: 

(i)  writing  down  the  physical  optics  integral  for  the  scattered  field 
without  any  shadowing  included; 

_ 

This  mfvjlfied  method  is  aomeiiineB  calied  extended  phyeicwl  optics". 
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(ii)  asymptotically  evaluatiag  this  e3Q)res8ioii  appropriate  to  an 
observation  point  at  lai^e  distances  from  the  sheet;  and 

(iii)  this  expression  to  a  discrete  angular  spectrum  of  waves 
to  obtain  their  amplitudes.  Since  these  amplitudes  are  independent  of  the  field 
point  we  have,  in  effect,  produced  an  exact  evaluation  of  the  int^ral  valid  cer¬ 
tainly  in  the  half  space  above  the  sheet. 

Let  us  take  the  equation  of  the  perfectly  conducting  sheet  to  be 

y^acosKx  (A.  i) 

where  a  is  the  amplitude  of  the  corrugafions  and  27r/K*d  is  the  period.  If 
the  incident  field  is  E-polarized,  we  write 

_E  =■  ze  (A.  2) 

(cf.  Eq.  2.  20),  implying 

V/  •  „A.  -jk{xsin0-ycos0)  /* 

H  = -YCcos0>:+sin0y)e  •’  (A,  3) 

where  Y  is  the  intrinsic  admittance  of  free  space. 

By  virtue  of  the  periodicity  of  the  sheet  and,  hence,  of  the  problem  as 
a  whole,  the  scattered  field  can  be  expressed  as  a  discrete  spectrum  of  waves 
which  waves  are  certainly  outgoing  as  regards  y  >  a.  Thus,  ws  have 


=  z 


5; — 1  -jk(xsin0  +vcose  ) 

>  A  e 

>  TR 

m*-a) 


(A.  4) 


where 


ksin0  *mK+ksin0, 
m 

kcosO  =  vk*” -(mK+ksinO)^ 

wy  * 


(A.  5) 


and  thafc  branch  of  the  square  root  is  diosen  having  imaginary  part  non-positive. 
Application  of  the  physical  optics  approximation,  followed  by  steps  (i)  throu^ 
(iii)  above  now  gives 


.  .m _ k 

I  5  ‘ 

yk  -(mK+ksin0/ 


„  mKa  .  „  i  ,  ,  , 

COS0- - sinS^J  (fx  ] 
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where 


ti  ^  z.fkcos6+  \/k^-(mK+ksin0)^  )  . 
rn  \  ^ 


In  particular. 


and  as  a  — >  0  , 


Aq  *  -Jg(2akcos  0) , 


A  1 ,  A  — >  0  ,  m  ^  0  , 
0  m 


in  agreement  with  the  known  solution  for  a  fiat  sheet. 
If  the  incident  field  is  H-polarized,  we  take 


implying 


„i  A  -jk{x  sin  0  -  y  cos  0) 
H  *ze 


i  1  /  -jk(xsin0-ycos0) 

E  =*  —  (cos0x+sin0y)e  ^ 


and  expand  the  scattered  magnetic  field  as 


A 

*  z 


03 


m=-oo 


A' 

m 


-jk(x8in0  +ycos0  ) 

^  Vrt  TYt 


where  0  is  as  before.  On  evaluation  of  the  physical  optics  integral, 
m 

obtain 


(A.  6) 


(A.  7) 


(A.  8) 


(A.  9) 


(A.  10) 

we  now 
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A*  *A  .  (A.  11) 

m  m 

The  equivalence  of  the  results  for  the  two  polarizations  is  consistent  with  &e 
known  fact  that  the  lAiysical  optics  approximation  for  a  perfect  conductor  is 
inherently  polarization  insensitive. 

Results  computed  from  Eqs.  (A.  6)  and  (A.  11)  agree  witii  those  given 
in  the  appropriate  columns  of  Tables  III-l  throng  ni-4.  Tables  V-l  throng 
V-4,  and  obtained  from  a  numerical  evaluation  of  the  integrals  in  Eqs.  (4. 4) 
and  (6.2). 


APPENDIX  B 

An  Alternative  Integral  Eqaaticn  for  E  Polarization 


Another  integral  equation  for  the  surface  current  on  a  periodic 
conducting  surface  can  be  derived.  It  is  different  from  Eq.  (2. 11),  and  since 
it  is  a  Fredholm  equation  of  the  second  kind  vdth  a  weakly  singular  keimel,  it 
can  be  solved  using  an  iteradvs  sdieme  such  as  a  Neumann  series. 

Consider  the  perfectly  conducting  t^Undrical  surface  of  Fig.  B-1,  The 
equation  of  the  surface  is  y  =  fix),  where  may  or  may  not  be  a  periodic 
function  of  x  at  this  stage  of  tlie  discussion.  The  incident  piano  electromag¬ 
netic  wave  propagates  in  the  direction  of  the  unit  vector 

k  =  sin0x-cos0y  ,  (B.  1) 


and  the  incident  electric  field  is  parallel  to  the  generators  of  the  surface: 

i_  -jk(x6in0-ycos0) 

E  —  e 

The  time  dependence  factor  e^  is  omitted. 

The  scattered  field  is  (Stratton,  1941) 


(B.2) 


where 


(B.3) 

{B.4) 


F(x,  y,  z)  and  P’(x’,y',z  ;  are  >he  field  and  integration  points  respectively, 

I_  is  the  volume  current  density,  and  v  is  a  volume  containing  both  the  field 
point  and  the  scattering  surface.  The  vectors  r  and  r'  lying  in  the  plane  z  =  0 
are  shown  in  Fig.  B-1  and  are  given  by 


r  =  xx+yy. 


t  I  1  f  ^ 

r'  *x'x+y'y 


(B.5) 
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21.1 


If  K  is  the  surface  current  density,  then 


Kds  =  Idv  , 


(B.6) 


and  therefore 


E®(r+zz) 

’-"7.  ~ 


^0  r 


\  ^ -  K  ds 

Js  ^ 

ip  00  -jkR. 
S*  J  -00  ^ 


(B.7) 


where  S'  is  the  scattering  surface,  and  S  is  its  profile  in  the  z  0  plane. 


SiOO  -jhRj 


(2) 

dz'  «  jTTHjj  '(kR) 


(B.8) 


where 


R  =  V^x-x»)  +(y-y’)  . 


(B.9) 


Substitution  of  Eq,  (B.8)  into  Eq.  (B.7)  fields 


E^(r+zz)  =  E®(x) 
z  —  z 


Js 


jfW)K  (r')ds' 
U  z~ 


(B.  10) 


where 


ds'  =  v  l+f'(3^)  dx'  . 


(B.li) 


(2)  f  2 

The  kernel  h'  (kR)  y  l  +  f  (x)  is  nonsymmetric  and  weakly  singular  at  R  =  0, 

"  (2) 
due  to  the  logarithmic  singularity  of  (kR). 

Let  the  normal  n  to  the  surface  S'  at  the  fixed  point  P  be  as  shown 

in  Fig,  B-2;  differentiating  both  sides  of  Eq.  (B.  10)  with  respect  to  n,  as  P 

approaches  the  surface  S',  we  obtain 
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aE®(r) 
z  “ 


P  on  S 


(r')d8» 
u  z  ~ 


(B.  12) 


To  evaluate  the  left-hand  side,  observe  that  with  the  three  basis  vectors 
t,  n  and  z  of  Fig.  B-2,  Maxwell's  equation  V<E  =  with  E  *E  gives 

■“  lr~  ”  Z 

(H  =  0  on  S') 
n 


aE  aE 


(B.  13) 


and  therefore. 


H  *0  , 
z 


(B.  14) 


Since  K  =  n  X  H  and  E  =  +  E®,  Eq.  (B.  14)  becomes 

—  zzz 


IT  ®  • 


(B.  15) 


jk  (f  sln0  -  f  COS0)  e-J“*s‘»l*-ycoB0) 


(B.  16) 


and  making  use  of  the  relations 


Ox  \  ^  ^  ^ 

■r"  =  (^x)  •  n  =  X  •  n  , 
an 


9y  iw,  \  ^ 

■^  =  (Vy)-n«yn  , 


vB.  1?) 


and  Eq.  (B.  1),  Eq.  (B.  16)  gives 
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z  .1  Z'  A  “jk(xsin0-ycos0) 
-  *  -jkk-n  e 


(B.  18) 


Substitution  of  Eqs.  (B.  15)  and  (B.  18)  into  Eq.  (b.  12)  yields 


K  (r).-  p  (^.-g-jk(xsin0-ycos0)^ 


h' 

4  an  J 


5fW)K  (r')ds'  . 
0  z  "■ 


(B.  19) 


Since  the  integrand  in  Eq.  (B.  19)  is  singular,  when  P  *  P’,  particular 

care  must  be  taken  in  interchanging  the  differentiation  with  respect  to  n  and 

the  integration  with  respect  to  x*.  This  difficulty  was  previously  encountered 

by  Maue  (1949)  and  Riblet  (1952),  among  others,  and  can  be  overcome  in  the 

manner  that  they  adopted.  The  procedure  consists  of  interchanging  9/ an  and 

\  ,  but  in  replacing  K  (r)  on  the  left-hand  side  with  (r):  for  a  detailed 
JS’ 

justification  of  this  step,  see  Maue  (1949).  The  result  is 


0  ,  -jk(x8in0-ycos  0)  , 
—  k*ne  + 


4^ 


(B.20) 


It  now  remains  to  evaluate  the  partial  derivative  in  the  integrand,  viz. 

3  0  .  3R  ,,, 


„(2Wr, 

a(kR)  '  "l  ■ 


where 


(B.22) 
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BR 

9ri 


dx  dn  9y  9n 


x-x' 
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•n+ Y^yn - ^ 


=  -R*n,  (B.23) 


with 


(B.  24) 


which  is  a  unit  vector  directed  from  the  fixed  point  P  towards  the  integrand 
point  P’.  When  Eqs.  (B.  21),  (B.  22)  and  (B.  23)  are  inserted  Into  Eq.  (B.  20), 
we  have 


K  (r) 
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„  Z'  A  -jk(xsin0-ycos0)  , 
•2  —  k'ne  + 

^0 


ik 
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Js 


(kR)K  (r')dB'  (B.26) 
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which  is  a  weakly  singular  one-dimensional  Fredholm  integral  equation  of  the 
second  kind  for  the  unknown  surface  current  density  K^. 

As  a  simple  check,  let  us  take  f(x)  =  0.  Then  R  *  n  *  0  and 


K 

z 


cos  0  e 


-jkx  sin  0 


(B.26) 


which  is  the  expected  geometric  optics  result. 

Equation  (B.  25)  is  valid  for  any  cylindrical  surface  and  thus  can  be 
applied  to  the  case  of  a  periodic  surface. 


